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1.  INTRODUCTION 


1.1  PROJECT  SUMMARY 

This  report  presents  the  results  of  the  Phased  Array  Laser  System 
(PALS)  project.  During  the  PALS  project,  analyses  and  design  approaches 
were  developed  for  a  space-based  HF  phase-conjugate  laser  system.  The 
preliminary  configuration  reviewed  in  this  report  is  the  result  of  a  number 
of  modeling  and  analytical  efforts,  and  some  subscale  testing  and  data 
review,  aimed  at  providing  the  first  step  toward  an  eventual  system  concep¬ 
tual  design. 

It  is  anticipated  that  HF  phase-conjugation  technology  development, 
made  possible  by  this  project,  may  become  the  basis  for  an  SDI  SBL  system 
design.  Such  a  system  has  substantial  advantages  of  design  and  performance 
margin  over  more  conventional,  non-phase-conj ugate  systems.  In  a  phase- 
conjugate  system,  device  beam  quality  is  limited  only  by  the  beam  quality 
of  a  small  oscillator  which,  given  the  low  oscillator  power  requirement, 
can  be  near-di f f raction  limited.  Furthermore,  the  key  requirement  of  mul¬ 
tiple  device  coupling  requires  minimal  alignment  control  (%5  waves  of  angu¬ 
lar  alignment,  and  ^3  meters  of  piston  control).  Similarly,  output  device 
jitter  depends  only  on  the  jitter  of  the  small  oscillator,  with  high  energy 
laser  ( HE L )  device  jitter  corrected  to  >1  MHz.  These  performance  and 
design  advantages  provided  powerful  motivation  for  the  PALS  project,  which 
has  developed  initial  approaches  to  the  design  of  such  a  system. 


gure 


2.  PROJECT  SUMMARY 


2.1  SYSTEM  REQUIREMENTS 

The  PALS  mission  requirements,  most  of  which  are  classified  (Ref.  2.1), 
are  based  on  the  SDIO  SBL  1985  technical  guidelines  document.  More  spe¬ 
cific  requirements  were  based  on  the  Lockheed  SBL  Phase  I IB  mission  analy¬ 
sis  report,  as  a  typical  example  of  SBL  brightness  levels.  However,  the 
PALS  phase-conjugate  design  has  a  key  advantage  over  conventional  HEL  sys¬ 
tems  for  brightness  scalability:  In  a  phase-conjugate  design,  complexity 
is  not  a  strong  function  of  system  power  and  diameter.  This  feature  makes 
the  PALS  approach  intrisically  scalable  to  much  higher  brightness  levels. 

2.2  SYSTEM  DESIGN  APPROACH 

The  performance  benefits  provided  by  nonlinear  phase  conjugation  for 
an  SBL  system  design  can  be  maximized  by  changing  the  design  at  the  system 
level.  Ultimately,  this  technology  makes  possible  two  important  areas  of 
performance  improvement: 

0  Transfer  of  virtually  the  complete  HEL  beam  quality  and  jitter 
budgets  from  the  HEL  device  to  a  low  power  oscillator  subsystem 

o  Phase-locking  of  multiple  HEL  gain  generator  outputs. 

The  system  configuration  that  achieves  this  is  a  Bidirectional  Phase 
Conjugate  MOPA.  As  shown  in  the  functional  block  diagram  in  Figure  2-1, 
there  are  four  subsystems: 

o  Oscillator  subsystem  o  Power  amplifier  subsystem 

o  Beam  director  subsystem  o  Phase  conjugation  subsystem 

The  oscillator  subsystem  provides  a  reference  beam  which  will  define  the 
beam  quality  and  jitter  of  the  full  HEL  system  after  later  amplification 
and  conjugation.  The  output  reference  beam  from  this  subsystem  propagates 
to  a  beacon  mirror,  which  would  be  the  key  element  on  a  very  small,  sepa¬ 
rate  beacon  satellite  used  only  for  that  purpose.  The  beam  leaving  the 
beacon  diverges  to  fill  the  composite  primary  of  a  single-aperture  tele¬ 
scope,  which  then  focuses  the  beam  to  a  demagnified  size  on  a  telescope 
secondary  mirror.  The  beam  at  this  point  includes  OPD  information  such  as 


2.1  "Final  PALS  Mission  Requirements  and  Guidelines  Document  (U)," 
TRW  AP-0382 ,  19  Feb  1986  (Secret) . 
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Figure  2-1.  Phased  array  laser  system 


the  relative  piston  errors  of  the  composite  primary.  After  propagation 
through  the  next  subsystem,  the  power  amplifier  subsystem,  these  and  other 
accumulated  OPD  errors  and  jitter  will  be  conjugated  in  the  phase  conju¬ 
gation  subsystem. 

As  the  conjugate  beam  returns  through  the  HEL  system,  it  corrects  for 
the  original  OPDs  by  passing  back  through  them  with  their  phase  conjugate. 
As  it  leaves  the  power  amplifier  subsystem,  the  beam  passe;  through  beam 
separation  optics  which  direct  the  beam  through  a  somewhat  different 
optical  path,  so  that  the  final  output  beam  is  approximately  collimated 
(directed  to  a  target)  after  correcting  primary  mirror  piston  errors. 
Figure  2-1  is  a  block  diagram  which  shows  the  relationship  between  the 
different  subsystems. 


2.3  SUBSYSTEM  DESIGN  APPROACHES  AND  ANALYSIS 


2.3.1  Master  Oscillator  and  Power  Amplifier  (MOPA)  Subsystems 

The  MOPA  activity  was  principally  concerned  with  development  of  pre¬ 
liminary  constraints  to  guide  selection  of  the  number  and  size  of  HEL 
amplifiers.  This  activity  also  involved  modeling  to  define  the  oscillator 
size  required  to  suppress  ASE  and  to  provide  good  extraction  from  the  power 
ampl ifiers. 

Preliminary  selection  of  appropriate  design  ranges  was  made  by  modi¬ 
fying  the  LFCM  amplifier  model  to  include  HF  kinetics.  This  model  was  used 
to  look  at  both  extraction  efficiency  and  ASE  issues  as  a  function  of  the 
relative  sizes  of  the  oscillator  and  amplifier  subsystems.  The  ANOM  and 
CLAM  models  were  also  used  to  look  at  the  effects  of  SBS  Doppler  shifting 
on  extraction  efficiency;  the  results  of  this  effort  showed  that  Doppler 
shift  compensation  using  a  near-sonic  axial  xenon  flow  is  unnecessary  to 
get  good  extraction. 

At  the  conclusion  of  this  modeling  effort,  review  of  the  results 
showed  that  the  required  oscillator  and  amplifier  sizes  are  well  within  the 
capability  ranges  of  HF  gain  generators.  This  result  was  based  on  assess¬ 
ments  of  the  power  needed  to  meet  three  requirements: 

1.  Adequate  margin  over  the  SBS  threshold  at  the  SBS  subsystem 

2.  Good  overall  amplifier  extraction  efficiency 

3.  ASE  suppression. 

Modularity  of  the  amplifiers  makes  them  particularly  well-suited  to  any 
further  power  scaling  that  may  be  necessitated  by  increases  in  requirements 
on  future  systems.  Figure  2-2  shows  the  power  amplifier  subsystem. 

Another  activity  included  in  this  effort  was  development  of  beam  sepa¬ 
ration  approaches.  Beam  separation  techniques  using  a  polarization  separa¬ 
tion  grating  rhomb  and  optically  active  coatings  were  modeled,  and  prelimi¬ 
nary  design  approaches  look  effective. 


2.3.2  Phase-Conjugation  Subsystem 

Utilization  of  the  BOUNCE  SBS  phase  conjugation  model  was  particularly 
key  to  development  of  design  approaches  for  this  subsystem.  Using  this 
model,  the  SBS  threshold  was  derived,  and  the  result  was  similar  to 
although  somewhat  lower  than  that  measured  in  NRL  experiments.  The  derived 
value  is  well-suited  to  the  SBL  operating  point. 

One  important  feature  of  the  SBS  subsystem  is  the  design  approach  and 
approximate  size  and  complexity  required  for  the  flowing  SBS  cell.  Due  to 
the  similarity  between  design  requirement  for  a  pulsed  design  and  a  CW 
design  flowing  cell,  it  was  possible  to  use  conceptual  designs  available 
GFE  from  an  HF  conjugation  experiment  design  program  (CLPCT).  Figure  2-3 
shows  features  of  the  preliminary  design  concept. 
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Figure  2-3.  PALS  preliminary  design  concept. 


Additional  work  done  in  this  area  included  some  subscale  testing  on 
basic  SBS  issues,  and  beam  combination  modeling  done  with  the  BOUNCE  code. 
The  modeling  effort  yielded  an  initial  design  point  for  multiple  beam 
combination  geometry,  using  40  atm  xenon,  as  discussed  below. 

2.3.3  Beam  Director  Subsystem 

The  baseline  beam  director  design  developed  on  this  project  was  aimed 
at  meeting  the  two  most  basic  design  constraints: 

1.  Bidirectional  operation  to  permit  conjugation  of  multisegment 
piston  errors 

2.  Reasonable  field  of  view  to  allow  pointing  within  a  target  field 
without  moving  the  primary. 

Conjugation  of  the  beam  director  is  a  powerful  tool  to  increase  performance 
margins  and  reduce  the  otherwise  severe  segment-to-segment  alignment  con¬ 
straints;  primary  mirror  conjugation  yields  n,2  orders  of  magnitude  improve¬ 
ment  in  piston  error  and  angular  alignment  tolerances. 

Both  of  the  above  design  constraints  were  met  by  the  design  approach 
developed  on  PALS.  The  basic  scheme  was  to  use  a  primary  with  a  conic 
constant  near  zero,  to  make  it  possible  to  maintain  the  same  primary/ 


secondary  optical  alignment  as  the  pointing  angle  changed.  A  modification 
of  this  approach  was  also  developed,  in  which  a  single  secondary  mirror 
replaces  earlier  approaches  which  required  separate  secondaries  for  low  and 
high  power  beams  {Figure  2-4). 

PRIMARY  MIRROR 


RHOMB 

Figure  2-4.  Beam  director  concept  showing  secondary  and  relay  optics 


3.  MISSION  REQUIREMENTS,  GUIDELINES,  AND  MISSION  ANALYSIS 


3.1  STARTING  POINT 


In  order  to  determine  the  kind  of  mission  for  which  PALS  is  to  be 
used,  existing  space-based  laser  programs  were  examined  to  provide  a 
technical  starting  point  for  the  PALS  program.  One  was  chosen  that  seemed 
representative  of  most  such  programs,  namely  the  Lockheed  Space-Based  Laser 
(SBL)  program.  That  program  was  examining  the  design  issues  associated 
with  a  high  brightness  HF  laser  much  like  the  proposed  PALS  design. 

The  SBL  program  made  use  of  a  technical  guideline  document  (Ref.  3.1) 
and  the  decision  was  made  to  use  this  document  for  the  PALS  program  also. 
This  document  specifies  an  SBL  mission  in  terms  of  1-on-P  performance 
guidelines,  which  detail  number  and  type  of  RVs  and  launch  vehicles  killed, 
and  gives  a  range  histogram  for  the  M-on-N  target  scenario,  complete  with 
trajectories,  launch,  and  target  coordinate  information.  Also  included  are 
target  vulnerability  parameters  for  the  specified  launch  vehicles  and  RVs, 
given  in  terms  of  minimum  intensity  and  fluence  requirements,  and  damage 
spot  size  requirements. 


Under  phase  1 1 A  and  I IB  of  this  program  Lockheed  performed  a  detailed 
mission  analysis  of  the  specific  threat,  and  determined  a  minimum  overall 
system  brightness  level  required  to  perform  the  mission.  Although  this  is 
likely  to  be  a  constantly  evolving  requirement,  continually  subject  to 
change  as  program  requirements  become  better  defined,  their  value  was  taken 
as  an  initial  design  point  for  this  important  parameter. 

The  definition  for  overall  brightness  is  based  on  the  Air  Force 
version  of  the  equation  and  is  shown  below. 


Brightness 


PA  (1 

-  e2)exp-(2 fit)2 

\2P2 

(0.431  X/D)2f(e). 

(3-1) 


3.1  "Space-Based  Laser  Phase  II  (U),”  USAF,  February  19,  1985. 


where 


P  is  system  output  power 
e  is  beam  obscuration 

l  is  phase  error  induced  by  the  beam  director 
X  is  wavelength 

<7j  is  jitter  beam  spread,  usually  expressed  in  terms  of  X/D 
A  is  aperture  area 

f(e)  is  correction  in  beam  spread  due  to  obscuration  and  is  equal 
to  1/(1  +  1.7  e1 2 3 4 5 6 7-8 9 10 11) 

is  diffractive  beam  spread  =  0.431  X/D 
D  is  primary  aperture  diameter 

p  is  output  beam  quality  assuming  a  perfect  beam  director 

It  is  important  to  note  that  most  of  the  key  system  parameters  appear 
in  this  equation.  Therefore,  if  one  parameter  is  changed,  the  equation 
dictates  what  other  parameter(s)  must  be  changed,  and  by  how  much,  in  order 
to  keep  overall  brightness  constant.  The  effect  of  some  of  these 
parameters  will  be  shown  later  in  this  section. 

3.2  MISSION  REQUIREMENTS 

The  mission  requirements  are  listed  in  Table  3-1.  Most  are  classified 
and  can  be  found  in  Reference  2.1. 

TABLE  3-1.  PALS  Mission  Requirements. 


1.  Minimum  damage  spot  diameter 

2.  Minimum  edge  fluence  requirement 

3.  Minimum  edge  intensity  requirement 

4.  Minimum  overall  system  brightness 

5.  Total  firing  time  capability 

6.  Maximum  target  engagement  time 

7.  Maximum  time  delay  to  reach  battle-ready  condition 

8.  Maximum  time  delay  for  laser  beam  initiation 

9.  Minimum  time  between  refuelings 

10.  Shut-off,  restart  capability 

11.  Acquisition,  pointing,  tracking  system  compatibility 


The  mission  derived  guidelines  are  presented  in  Table  3-2.  Most  of 
these  are  internally  generated  system  parameters  which  are  intended  to 
guide  the  design,  and  may  be  changed  more  frequently  than  the  requirements, 
if  this  becomes  necessary. 

TABLE  3-2.  PALS  Mission  Derived  Guidelines. 

1.  Minimum  system  output  power 

2.  Primary  mirror  diameter 

3.  Maximum  overall  system  jitter 

4.  Maximum  allowable  beam  director  large  angle  slew  acceleration 

5.  Minimum  beam  director  large  slew  angle  capability 

6.  Minimum  retargeting  time  within  a  target  field 

7.  System  beam  quality  assuming  a  perfect  beam  director,  stated 
as  a  range  between  a  minimum  and  maximum  value 

8.  Beam  director  wavefront  error 

9.  Weight  and  volume  optimization  as  system  consideration 

10.  Design  for  operational  and  nonoperational  conditions 

11.  Minimum  overal 1  lifetime 

12.  Survivability  in  specified  battle  conditions 

13.  Simplicity  and  reduced  complexity 


3.3  MISSION  ANALYSIS 

In  order  to  determine  the  effect  of  the  major  terms  in  the  brightness 
equation,  Equation  (3-1)  is  rewritten  as  follows: 

Brightness  =  —  ^  .  exp-(2*^)2  .  ~  .  - K, - 

*  P  o} 

1  +  - J - p - 

(.431  X/Dr  f(e) 
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It  is  of  interest  to  determine  the  effect  of  the  three  coefficients  of 
the  PA/>2  term.  This  is  shown  in  Figure  3-1.  In  order  to  keep  the  overall 
brightness  greater  than  the  required  value,  the  product  of  the  three 
coefficients  must  be  greater  than  0.61.  That  is  a  decrease  in  one  factor 
must  be  compensated  by  an  increase  in  one  or  both  of  the  other  terms. 

The  above  values  are  for  optimum  range  conditions.  At  ranges  greater 
than  the  optimum,  the  intensity  in  the  damage  spot  decreases  with 
increasing  range  due  to  diffractive  effects.  Figure  3-2  shows  this  effect 
for  peak,  average,  and  edge  intensities  over  the  ranges  of  interest.  All 
values  have  been  normalized.  Note  that  at  ranges  less  than  the  optimum 
range  the  intensities  remain  constant,  due  to  the  defocus  of  the  beam  nec¬ 
essary  to  provide  the  required  intensity  at  the  edge  of  the  specified 
damage  spot  (otherwise  the  spot  would  be  smaller  than  the  required  damage 
spot  size).  An  alternate  method  for  accomplishing  this  is  to  scan  the  beam 
around  the  edge  of  the  spot,  however,  this  results  in  a  decrease  in  inten¬ 
sity  in  the  center  of  the  spot,  which  might  be  undesirable. 

When  examining  the  two  vulnerability  requirements,  edge  intensity  and 
fluence,  it  is  useful  to  determine  which  one  has  the  dominant  effect  as 
ranges  are  increased.  The  answer  is  that  the  edge  fluence  is  the  more 
stringent  of  the  two,  and  the  target  dwell  time  must  be  increased  to  keep 
the  edge  fluence  at  or  above  the  minimum  required  level.  Figure  3-3  shows 
the  effect  on  target  dwell  time  over  the  ranges  of  interest.  All  values 
have  been  normalized. 


0  - 1 - 1 - 1 - 1 - 1 - T - T - T - 1 - 1 - 1 - 1  I  ! 

0.1  0.3  OS  0.7  0.9  1.1  1.3  IS 

- —  OR  2tt0  OR  P 
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Figure  3-1.  Relative  brightness  affected  by  .jitter  WFE,  BQ,  and  OBSC 
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Figure  3-3.  Edge  intensity  and  dwell  ti 


4.  PALS  SUBSYSTEM  OVERVIEW 


4.1  FUNCTIONAL  DESCRIPTION 

The  Phased  Array  Laser  System  (PALS)  preliminary  configuration  has 
been  chosen  to  satisfy  the  mission  requirements  and  guidelines  (classi¬ 
fied),  as  well  as  the  system  and  subsystem  requirements  listed  in  the  PALS 
project  requirements  document  (PRD).  The  preliminary  configuration  has 
been  divided  into  four  subsystems  and  nine  assemblies  listed  below: 

1.  Master  Oscillator  Subsystem 
Oscillator  assembly 

Oscillator  output  optical  train  assembly 

2.  Power  Amplifier  Subsystem 
Ampl i f ier  assembly 

Beam  separation  and  recombination  assembly 
Amplifier  optical  train  assembly 

3.  Beam  Director  Subsystem 

Beam  director  assembly 
Beacon  assembly 

4.  Phase-Conjugation  Subsystem 

Conjugation  optical  train  assembly 
Conjugation  assembly 

These  subsystems  are  shown  schematically  in  Figure  4-1,  and  function¬ 
ally  in  Figure  4-2.  In  this  configuration  the  beam  begins  at  the  oscil¬ 
lator  subsystem  at  the  left  of  Figure  4-1  and  is  directed  to  the  free- 
fiying  beacon  mirror  which  is  off  the  right  side  of  the  page  at  a  distance 
of  approximately  20  times  the  primary  mirror  diameter  away  from  the  primary 
mirror.  The  beam  is  then  directed  back  towards  the  primary  mirror, 
diverging  to  fill  the  entire  aperture. 

The  low  power,  near  dif fraction-1 imi ted  beam  is  then  directed  towards 
the  secondary  mirror  which  sends  it  towards  the  amplifiers  via  the  beam 
director  pointing  optics  and  amplifier  separation  and  recombination  optics. 
The  smaller  beams  are  directed  to  the  individual  amplifiers,  where  they  are 
amplified  and  directed  to  the  phase-conjugation  subsystem. 
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The  phase-conjugation  optical  train  spectrally  separates  the  beams 
into  separate  lines  and  directs  them  to  the  SBS  cells.  A  conjugate  beam  is 
reflected  back  from  the  cells  and  propagates  back  along  the  input  beam 
path.  It  is  spectrally  recombined  and  sent  to  the  amplifiers  for  amplifi¬ 
cation  to  maximum  power. 

The  high-power  beam  is  then  directed  towards  the  beam  director 
pointing  optics  after  passing  through  the  isolation  grating  rhomb  which 
spatially  separates  the  high-  and  low-power  beams  due  to  their  different 
polarizations.  The  high-power  beam  is  then  directed  to  the  secondary  and 
primary  mirrors  from  which  it  is  propagated  to  the  target. 

4.2  SUBSYSTEM  DESCRIPTION 

4.2.1  Oscillator  Subsystem 

The  primary  function  of  the  oscillator  subsystem  is  to  provide  a  very 
good  beam  quality,  linearily  polarized  low-power  reference  beam  to  the 
amplifiers  (via  the  beacon  and  primary  mirror). 

This  is  accomplished  by  using  a  small  HF  laser  with  a  cylindrical  gain 
generator  much  like  the  larger  amplifier  design.  The  oscillator  beam  must 
be  spectrally  matched  to  the  power  amplifiers  to  provide  good  extraction 
efficiency  and  to  limit  to  four  the  number  of  lines  on  which  the  major 
portion  of  the  output  power  is  contained.  It  also  must  start  with  suffi¬ 
cient  power  to  propagate  the  relatively  long  distance  to  the  beacon  mirror 
and  still  provide  the  required  input  power  to  the  amplifiers. 

The  oscillator  optical  train  will  include  optics  to  remove  high  fre¬ 
quency  jitter  from  the  beam  and  direct  it  to  the  beacon  within  the  allow¬ 
able  pointing  requirements.  An  enlargement  of  the  oscillator  subsystem  is 
shown  in  Figure  4-3. 

4.2.2  Power  Amplifier  Subsystem 

The  primary  function  of  the  power  amplifier  subsystem  is  to  receive 
the  oscillator  beam,  amplify  it  bidirectionally,  and  delivery  the  high 
power  beam  to  the  beam  director. 

This  is  accomplished  by  using  seven  scaled-up  ALPHA  II  type  amplifiers 
arranged  in  a  parallel  configuration  with  the  exhaust  gases  directed 
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Figure  4-3.  Cylindrical  oscillator. 


radially  outward  to  avoid  contamination  of  optics  or  absorption  of  the 
beam.  The  amplifiers  are  operated  bidirectionally  in  order  to  phase  conju¬ 
gate  most  of  the  gain  generator  aberrations.  Beam  separation  and  recombi¬ 
nation  optics  are  provided  to  direct  the  beam  to  and  from  the  amplifiers, 
with  coarse  correction  of  low-power  beam  alignment  errors.  An  isolation 
grating  rhomb  is  used  to  spatially  separate  the  high-  and  low-power  beams 
due  to  the  difference  in  polarization  caused  by  the  quarter-wave  coating  in 
the  phase  conjugation  subsystem.  An  enlargement  of  the  power  amplifier 
subsystem  is  shown  in  Figure  4-4. 

4.2.3  Beam  Director  Subsystem 

The  primary  function  of  the  beam  director  subsystem  is  to  receive  the 
low-power  oscillator  beam  and  direct  it  to  the  amplifiers,  then  receive  the 
amplified  high-power  beam  and  direct  it  to  the  target,  with  the  necessary 
pointing  accuracy. 

This  is  accomplished  using  a  primary,  secondary,  tertiary,  and  relay 
mirrors,  as  well  as  a  free-flying  lightweight  beacon  mirror  to  relay  the 
low  power  oscillator  beam  to  the  amplifiers.  The  primary  mirror  is  spher¬ 
ical,  gimbaled,  mul ti segmented  (with  identical  segments),  and  has  piston 
and  jitter  errors  corrected  by  the  phase-conjugation  subsystem.  A  single 
secondary  is  used  for  both  the  high-  and  low-power  beams  since  this  mini¬ 
mizes  obscuration  and  enables  the  secondary  to  be  conjugated.  Small  angle 
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The  primary  function  of  the  phase-conjugation  subsystem  is  to  accept 
the  partially  amplified  beam  from  the  amplifiers,  and  reflect  a  phase-con¬ 
jugated  beam  back  through  the  amplifiers. 

This  is  accomplished  by  using  a  flowing  SBS  cell  concept  similar  in 
design  to  the  cell  to  be  used  in  the  COS  experiment  of  the  CLPCT  program. 
This  cell  utilizes  xenon  at  40  atmospheres  pressure  and  300  K,  which  is 
filtered,  purified  and  flowed  transversely  across  the  beam  at  between  20 
and  80  m's.  Two  helium  aerowindows,  one  each  at  the  entrance  and  exit 
apertures,  are  used  to  confine  the  xenon  within  the  SBS  cell  at  the  correct 
pressure,  with  minimum  beam  distortion.  Individual  SBS  cells  are  used  for 


Figure  4-5a.  Beam  director  concept  showing  secondary  and  relay  optics. 
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Figure  4 - 5b .  Free-flying  beacon  mi rror--re1 ays  oscillato 
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Figure  4-6.  Functional  diagram  of  beam  director  subsystem. 

each  spectral  line,  since  it  will  probably  be  simpler  to  construct  several 
small  cells  than  one  large  one,  and  also  to  prevent  the  focal  regions  from 
disturbing  each  other.  Note  that  the  focal  volume  of  each  cell  contains 
the  combined  energy  from  all  the  amplifiers  contained  on  a  single  line.  An 
enlargement  of  the  SBS  cell  configuration  showing  four  side-by-side  cells 
is  shown  in  Figure  4-7. 

4.3  PRELIMINARY  POWER  BUDGET 

Table  4-1  is  a  preliminary  power  budget  used  to  determine  the  oscil¬ 
lator  power  required  to  supply  the  output  power  specified  in  reference  4.1. 
Best  possible  estimates  were  used  for  mirror  reflectivities  and  grating 
efficiencies,  and  propagation  loss  to  the  beacon  was  based  on  an  initial 
beam  quality  of  1.2  and  a  bucket  size  of  0  =  0.8.  Amplifier  gain  was  taken 
as  the  ratio  of  the  output  power  to  the  power  reflected  from  the  SBS  cell, 
which  is  also  equal  to  the  ratio  of  the  SBS  input  power  to  the  input  power 
supplied  by  the  oscillator.  Both  were  calculated  for  ideal  cases  assuming 
no  losses  in  the  system. 


4.1  Final  PALS  Mission  Requirements  and  Guidelines  Document  (U),  AP-0382, 
February  19,  1986  (Secret). 


21 


XENON  EXHAUST 


NOZZLE 


MIRROR 


HELIUM 
EXHAUST  A 
OEFLECTOR 


FOCAL 

REGION 


XENON  TANK 


HELIUM  TANK 


INPUT/OUTPUT 
TO  AMPLIFIERS 


Side  View  Top  View 

Figure  4-7.  Transverse  flowing  SBS  cell  concept. 
TABLE  4-1.  Preliminary  PALS  power  Budget. 


Forward 

Pass 


Component 

Transmission 

Normalized 

Power 

Oscillator 

0.150 

Isolation  optic 

0.995 

0.149 

Propagation  loss 

0.550 

0.082 

Beacon  relay 

0.985 

0.081 

Primary 

0.995 

0.080 

Secondary 

0.995 

0.080 

Six  relay  optics 

0.970 

0.078 

Isolation  rhomb 

0.640 

0.050 

Beam  distribution 

0.995 

0.049 

Amplifier 

32.300 

1.598 

SBS  grating 

0.980 

1.566 

Quarter  wave 

0.970 

1.519 

Focusing  mirror 

0.995 

1.511 

SBS  cell 

0.680 

1.028 

Focusing  mirror 

0.995 

1.023 

Quarter  wave 

0.970 

0.992 

SBS  grating 

0.980 

0.972 

Amplifier 

32.300 

31.397 

Beam  distribution 

0.995 

31.240 

Isolation  rhomb 

0.980 

30.615 

Six  relay  optics 

0.970 

29.708 

Secondary 

0.995 

29.560 

EACH  LINE 
FOCUSES  TO  A 
DIFFERENT  CELL 


5.  MOPA  (MASTER  OSCILLATOR  AND  POWER  AMPLIFIER)  SUBSYSTEMS 


5.1  MOPA  SUBSYSTEM  PRELIMINARY  CONCEPT  AND  TRADE  STUDIES 

The  objective  of  this  task  was  to  develop  analytical  tools  to  optimize 
the  MOPA  configuration  with  respect  to  mission  requirements  and  guidelines. 
The  approach  was  to  perform  parametric  studies  on  ALPHA-type  cylindrical 
amplifiers  to  determine  a  preliminary  design  point.  The  parameters  to  be 
optimized  include  the  number  of  amplifiers,  the  amplifier  length,  diameter, 
input  power,  output  power,  and  power  to  the  SBS  cell.  Issues  to  be  con¬ 
sidered  include  suppression  of  amplified  spontaneous  emission  (ASE),  and 
the  interaction  of  gain  saturation  in  the  amplifier  with  phase  conjugation, 
which  depends  on  the  amplifier  Fresnel  number  and  the  saturation  level. 

Several  analytical  models  have  been  used  to  address  these  issues. 

They  include  the  ANOM  and  CLAM  codes,  which  are  one-dimensional  models  of 
CW  HF  resonators  and  amplifiers,  respectively,  and  are  anchored  to  the 
three-dimensional  model  of  the  ALPHA  device  (the  CROQ  code).  Also  used  was 
the  LFCM  code,  which  is  a  two-dimensional  model  of  an  amplifier  and 
includes  an  ASE  model.  These  codes  are  described  in  more  detail  in  the 
next  section. 

The  preliminary  configuration  for  the  gain  generator  geometry  was 
chosen  to  provide  the  required  power  with  the  minimum  design  risk.  LFCM 
was  used  to  perform  a  systems  trade  between  the  gain  generator  diameter, 
the  gain  length,  and  the  number  of  amplifiers.  The  constraints  were  that 
the  total  output  power  must  equal  the  PALS  requirement,  and  that  the  gain 
generator  diameter  and  length  must  be  less  than  the  structural  and  align¬ 
ment  limits  set  by  the  ALPHA  program.  The  normalized  input  power  was 
0.1/(number  of  amplifiers)  kW,  and  the  SBS  reflectivity  r,as  0.68.  The 
results  of  these  calculations  are  shown  in  Figure  5-1.  An  additional 
consideration  in  the  choice  of  configuration  is  that  while  fewer  amplifiers 
minimize  the  system  complexity,  a  greater  gain  generator  volume  is  required 
for  each  amplifier.  The  preliminary  configuration  choice  is  shown  in 
Figure  5-1  and  consists  of  seven  amplifiers  9  m  long  and  2.7  m  in  diameter. 
The  arrangement  of  these  seven  amplifiers  is  shown  in  section  view  in 
Figure  5-2  and  in  lateral  view  in  Figure  5-3. 
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Figure  5-1.  Gain  generator  geometry  chosen  to  provide  required  power 
with  minimum  design  risk. 
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Figure  5-3.  MOPA  subsystem  concepts. 
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More  detailed  analyses  on  components  of  the  MOPA  subsystem  are 
described  in  the  following  sections.  These  include  the  MOPA  extraction 
efficiency  and  beam  quality  and  beam  separation  issues  such  as  design  of 
the  grating  rhomb  and  the  quarter-wave  retardation  coatings. 

5.2  MOPA  EXTRACTION  EFFICIENCY  AND  BEAM  QUALITY* 

The  overall  objective  of  the  MOPA  analysis  task  was  to  develop  the 
analytical  tools  needed  to  define  the  MOPA  configuration.  In  particular, 
these  tools  will  be  used  to  define  the  conceptual  changes  to  modify  the 
ALPHA  resonator  design  to  an  amplifier  design.  The  analytical  models 
developed  are  intended  to  address  the  following  issues: 

•  Master  oscillator  power  and  bandwidth  required  to  provide  optimum 
extraction  efficiency,  good  beam  quality,  amplified  spontaneous 
emission  (ASE)  suppression,  and  adequate  SBS  reflectivity  in  the 
power  amplifier. 

•  Effect  of  the  SBS  frequency  shift  on  amplifier  extraction 
efficiency. 

•  Effect  of  amplifier  gain  saturation  on  the  ability  of  phase 
conjugation  to  correct  for  aberrations. 

The  approach  to  this  task  was  to  modify  one-  and  two-dimensional 
models  of  CW  HF  lasers.  The  ANOM  and  CLAM  codes  are  one-dimensional  models 
of  a  resonator  and  an  amplifier,  respectively.  They  both  contain  the  ACLOS 
gain  model,  which  describes  the  effect  of  both  inhomogeneous  and  homo¬ 
geneous  broadening  in  cascade-connected  transitions  of  a  CW  HF  laser.  This 
gain  model  is  based  on  pumping  and  quenching  rates  anchored  to  a  full 
kinetics  model  of  the  ALPHA  gain  medium  including  the  effects  of  rotationa1 
nonequilibrium  and  mixing  in  a  flowing  gain  medium.  The  ALPHA  gain  model 
has  been  validated  by  comparison  with  data  from  the  ALPHA  verification 
module  on  small  signal  and  saturated  gain,  output  and  intracavity  power 
spectra,  and  output  intensity  profiles  (reference  5.1).  The  output  of  the 
ANOM  and  CLAM  codes  includes  intensity  on  each  longitudinal  mode  and 
saturated  gain  and  index  as  a  function  of  frequency. 

5.1  SCORE  Review  Report  (U),  TRW  35980-6077-SX-00,  22  July  1983  (Secret). 
*This  section  reports  work  partly  done  on  the  Coupled  Amplifiers  program. 


The  LFCM  code  is  a  two-dimensional  model  of  an  amplifier  which  can  be 
used  to  predict  amplifier  beam  quality  as  well  as  extraction  efficiency. 

LFCM  was  originally  developed  for  excimer  systems.  Linder  PALS  it  has  been 
modified  to  model  the  ALPHA  amplifier  gain  generator  configuration  and  the 
ALPHA  gain  distribution  in  the  flow  coordinate.  Modifications  have  been 
made  to  permit  input  of  phase  screens  describing  the  medium  aberrations  as 
calculated  by  fluid  dynamics  codes.  LFCM  can  also  accept  input  fields  from 
the  wave  optics  model  of  the  optical  train  which  include  the  effects  of 
primary  mirror  segment  piston  and  tilt. 

5.2.1  ALPHA  Device  Amplifier  Models  in  LFCM 

Three  configurations  have  been  implemented  in  the  LFCM  (large  Fresnel 
number  conjugated  MOPA)  amplifier  code  to  model  extraction  efficiency  and 
beam  quality  in  the  amplifier  subsystem.  The  linear  geometry  is  redescribed, 
the  newly  added  "two-pass"  geometry  is  explained  and  input  parameters  of 
interest  are  defined  (although  explanations  of  these  parameters  are  also 
found  in  the  comment  section  at  the  beginning  of  the  source  code  and  in  ref¬ 
erence  5.2). 

5.2. 1.1  Cavity  and  Gain  Length  Geometry 

Figure  5-4  shows  a  schematic  of  an  ALPHA  resonator.  In  this  ring 
resonator,  the  circulating  beam  enters  at  the  waxicon,  is  annularized  to 
the  waxicon  outer  cone,  expands  as  it  travels  to  the  rear  element,  becomes 
smaller  as  it  goes  to  the  reflaxicon  outer  cone,  and  is  recompacted;  the 
beam  then  gets  scraped  and  the  feedback  portion  travels  to  the  waxicon 
inner  cone.  The  beam  passes  through  the  annular  leg  (and  through  the  gain 
medium)  twice. 

The  ALPHA  resonator  cofiguration  has  been  modified  slightly  to  produce 
a  "two-pass"  amplifier.  The  beam  is  assumed  to  be  annularly  symmetric,  so 
that  LFCM  can  be  reasonably  accurate  even  though  it  models  only  one  annular 
|  slice. 

!  LFCM  models  the  two-pass  amplifier  shown  in  Figure  5-5,  when  the 

•  logical  variable  ALPHA2  is  set  to  TRUE.  Figure  5-5a  shows  LFCM's  "cavity" 


5.2  "Geometry  Conf i gurations  for  the  MOPA  Code,"  D.  Copeland  and 
C.C.  Shih,  TRW  83 . K323-43 ,  1983. 
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Figure  5-5b.  Gain  and  field  sheets  in  the  LFCM  code. 

and  relates  it  to  the  ALPHA  device.  Figure  5-5b  shows  the  positions  of  the 
field  and  gain  sheets  and  their  numbering.  The  inputs  to  LFCM  are  L,  the 
cavity  length;  LG,  the  gain  length  (assumed  equal  to  L  unless  input);  and 
GAP1,  the  distance  from  the  outer  cones  to  the  beginning  of  the  gain 
medium.  The  gain  generator  is  assumed  to  be  centered  between  the  rear 
element  and  the  outer  cones. 

If  Rqg  is  the  radius  of  the  gain  generator  (cm),  then  L  =  Rqq  +  3  cm  + 
GAP1  -  LG  +  GAP  1 .  A  nominal  value  of  40  cm  for  GAP]  has  been  used  in  the 
calculations. 

When  the  gain  medium  is  shorter  than  the  cavity  (LG  <  L),  there  are 
some  gain  sheets  where  the  gain  medium  does  not  fill  the  space  between  the 
field  sheets  or  is  absent  altogether.  In  order  to  handle  this,  the  code 
defines  a  gain  length  window  function,  WZ,  which  consists  of  NZ  numbers, 
each  of  which  is  the  fraction  of  space  between  corresponding  field  sheets 
filled  by  gain  medium.  (For  the  case  illustrated  in  Figure  5- 5b ,  the 
values  of  WZ  would  be  0,  .65,  1,  1,  ...,  1,  .4,  0.)  The  gain  is  calculated 
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at  each  gain  sheet,  then  multiplied  by  the  window  function,  so  that  the 
appropriate  amount  of  gain  is  applied  to  the  field. 

In  this  "two-pass"  amplifier  configuration,  there  are  three  options. 
The  first,  with  the  parameter  ONEWAY  true,  sends  the  beam  into  one  axicon, 
through  the  gain  medium  twice,  and  then  out  the  other  axicon.  There  is  no 
phase  conjugator  or  mirror.  In  the  second  and  third  options,  with  ONEWAY 
false,  the  beam  is  sent  into  one  axicon,  through  the  gain  medium  twice, 
then  out  the  other  axicon;  the  beam  is  reflected  by  a  mirror  (MIRROR  true) 
or  by  a  perfect  phase  conjugator  (PERPHS  true),  returns  twice  again  through 
the  gain  medium  to  the  axicon  by  which  it  originally  entered,  and  then  is 
output.  (With  MIRROR  true,  the  reflectivity  is  100%;  with  PERPHS  true,  the 
amplitude  reflectivity  RC0NJ6  is  input.)  The  first  and  second  options 
enable  modeling  of  an  amplifier  in  a  system  using  "conventional"  optics; 
the  third,  "nonlinear"  optics.  The  first  option  also  allows  for  cali¬ 
bration  of  the  amplifier  model  by  simulation  of  one  round  trip  through  an 
oscillator  and  then  comparison  of  the  results  to  those  of  the  wel 1 -anchored 
oscillator  models  (such  as  BLAZER  and  CROQ) . 

The  choice  of  whether  the  beam  enters  the  annular  leg  at  the  waxicon 
or  the  ref  1 axicon  is  controlled  by  the  parameter  FORWARD.  Since  in  the 
ALPHA  resonator  the  beam  enters  at  the  waxicon,  this  will  correspond  to 
FORWARD  true.  FORWARD  false  will  tell  LFCM  to  have  the  beam  enter  at  the 
ref 1 axi con. 

In  the  case  of  a  "linear"  amplifier  (LINEAR  true),  the  device  is  as 
shown  in  Figure  5-6a.  The  beam  enters  at  one  axicon,  a  collimated  beam 
passes  through  the  annular  leg,  then  the  beam  exits  at  the  other  axicon. 

The  ONEWAY,  MIRROR,  and  PERPHS  options  work  for  this  configuration  also,  as 
well  as  for  the  previously  used  multiple-round-trip  (NRT  >  1)  options. 

5.2. 1.2  Transverse  Geometry  --  Linear  Amplifier 

Figure  5-6b  shows  the  geometry  of  the  linear  amplifier  as  seen  by  LFCM. 
WX  is  the  total  field  width,  including  guard  bands.  WX  =  OVERSIZ  *  DR.  DR 
is  the  geometric  size  of  the  right-hand  mirror.  The  extent  of  the  gain 
medium  (GX)  is  equal  to  WX,  unless  otherwise  input.  DL  is  the  size  of  the 
left-hand  (beam  input)  mirror;  it  equals  DR  unless  otherwise  input.  In 


LFCM,  the  beam  always  enters  on  the  left.  For  this  linear  case,  the 
transverse  coordinate  X  will  have  the  value  zero  in  the  center  of  the 
transverse  direction,  as  shown  by  the  dot-dash  line  in  the  figure. 

There  are  transverse  window  functions  which  make  sure  the  fields  and 
the  gain  roll  off  smoothly  but  quickly  to  zero  outside  of  their  geometric 
apertures.  These  are  set  by  calls  in  the  code  to  subroutine  WINDOW. 


The  amplified  spontaneous  emission  (ASE)  model,  described  in  reference 

5.3  by  Shih,  has  been  modified  to  work  for  ONEWAY  true. 

5. 2. 1.3  Transverse  Geometry  --  ALPHA  “Two-Pass'*  Amplifier 

Figures  5-7a  and  5- 7b  show  the  geometry  of  the  ALPHA  device  two-pass 
amplifier  model.  (FORWARD  false  has  been  assumed  in  the  figures.)  The  new 
parameters  are: 

OFFSET  Distance  between  the  gain  generator's  nozzle  exit  plane  and 
the  inner  geometric  edge  of  the  field.  Default  =  0.3  cm. 
LFCM's  X=0  is  at  the  nozzle  exit  plane  (NEP) . 

WXOCONE  Waxicon  outer  cone  radius.  Default  =  1.0  cm. 

RFOCONE  Reflaxicon  outer  cone  radius.  Default  =  1.3  cm. 

Other  parameters  are: 

RCONE  Rear  element  radius.  Default  =  2.6  cm. 

WX  Largest  field  size;  WX  is  set  to  RCONE  *  OVERSIZ.  (OVERSIZ 

must  be  large  enough;  about  1.8  is  recommended.) 

FX  Input  field  size;  FX  is  set  to  (WXOCONE  or  RFOCONE)  * 

OVERSIZ. 

GX  The  extent  of  the  gain  region;  GX  is  set  to  RCONE  +  OFFSET. 

5.2.2  Master  Oscillator  Power  Requirements 

Calculations  performed  on  PALS,  with  support  from  results  of  the 
Coupled  Amplifiers  program,  were  used  to  set  the  requirements  on  master 
oscillator  power.  These  calculations,  performed  with  the  ANOM  and  CLAM 
models,  demonstrated  that  the  power  requirement  at  the  phase  conjugator 
could  be  met  simultaneously  with  meeting  a  requirement  that  the  power 
extracted  from  a  bidi rectional  amplifier  be  approximately  the  same  as  that 
extracted  from  a  resonator  with  the  same  size  gain  generator.  These 
calculations  included  the  effect  of  the  SBS  frequency  shift  (see  below). 


5.3  "ASE  Model  in  MOPA  Code,"  C.C.  Shih,  TRW  83.K323-69  (and  83.K326-73 
for  details  of  derivation  of  F)  (FA  in  the  code),  1983. 
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LFCM  for  ALPHA2=.TRUE.  and  forward  =.FALS 


The  objective  of  this  task  was  to  determine  if  the  SBS  frequency  shift 
significantly  affected  bidirectional  amplifier  extraction  efficiency  caused 
by  the  offset  of  the  return  beam  from  the  line  center,  and  whether  flowing 
the  SBS  medium  along  the  direction  of  Stokes  beam  propagation  to  decrease 
the  SBS  shift  would  be  necessary  to  improve  the  extraction  efficiency.  The 
analysis  performed  with  the  ANOM  and  CLAM  codes  demonstrates  that 
decreasing  the  SBS  frequency  shift  is  not  needed  for  xenon  SBS  to  obtain 
the  same  extraction  efficiency  from  a  bidirectional  MOPA  as  from  a 
resonator. 

The  SBS  frequency  shift  obtained  from  40  atmospheres  of  xenon  is  107 
MHz,  while  the  unsaturated  g.'in  has  a  full  width  at  1/e  of  480  MHz  (a  unit 
Doppler  width  equals  240  MHz  in  the  figures).  The  input  to  the  amplifier 
model  (CLAM  code)  is  the  output  1 ongi tudi nal -mode  power  spectrum  from  the 
resonator  model  for  the  master  oscillator  (ANOM  code).  The  saturation 
level  in  the  master  oscillator  has  been  adjusted  so  that  the  output 
bandwidth  of  the  longitudinal  mode  spectrum  provides  optimal  saturation  and 
extraction  of  the  gain  in  the  power  amplifier.  Amplifier  extraction  was 
calculated  for  the  case  with  the  full  SBS  frequency  shift  (identified  as 
shifted),  and  for  the  case  where  the  SBS  frequency  shift  was  completely 
cancelled  by  Mach  1  flow  of  the  SBS  medium  in  the  direction  of  Stokes 
propagation  (unshifted).  The  output  intensities  for  these  two  cases  for 
two  cascade  connected  lines  are  compared  in  Table  5-1.  The  output  inten¬ 
sity  as  a  function  of  longitudinal  mode  is  shown  in  Figure  5-8,  and  the 
saturated  gain  is  plotted  in  Figure  5-9.  Since  the  oscillator  can  be 
designed  to  give  a  broad  spectrum  of  longitudinal  modes  above  threshold, 
the  amplifier  gain  is  wel:  saturated  in  the  wings,  and  parasitics  and  ASE 
should  not  be  a  problem. 

5.2.4  Effect  of  Gain  Distribution 

The  small  signal  gain  distribution  g0  as  a  function  of  the  flow  coor¬ 
dinate,  x,  for  the  P?(7)  line  in  the  ALPHA  device  is  shown  in  Figure  5-10. 
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TABLE  5-1.  High-Speed  Flow  of  SBS  Xenon  Medium  Unnecessary. 
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The  zero  crossing  occurs  at  ~2.3  cm  from  the  nozzle  exit  plane.  In  the 
current  ALPHA  design,  the  beam  in  the  annular  leg  expands  from  1.3  cm  at 
the  entrance  to  2.6  cm  at  the  rear  cone  with  a  0.3  cm  offset  from  the 
nozzle  exit  plane.  Consequently,  the  gain  profile  significantly  affects 
the  output  intensity  distribution,  producing  the  narrowed  profile  shown  in 
Figure  5-11.  The  calculation  was  done  for  propagation  through  the  annular 
leg  only,  assuming  a  uniform  intensity  beam  as  input  to  the  annular  leg. 

In  the  ALPHA  amplifier  configuration,  the  use  of  axicons  to  compact  and 
expand  the  beam  between  the  compact  and  annular  legs  produces  a  nonuniform 
intensity  distribution  as  input  to  the  annular  leg.  The  appropriate 
choices  of  waxicons  and/or  reflaxicons  and  optimization  of  the  beam  widths 
in  the  annular  leg  should  result  in  a  reasonably  uniform  intensity  distri¬ 
bution  in  the  compact  leg. 

5.2.5  Interaction  of  Gain  Saturation  with  Phase  Conjugation 

LFCM  was  used  to  determine  the  magnitude  of  the  interaction  between 
gain  saturation  and  phase  conjugation  in  amplifiers  of  7,  9,  and  11  m  gain 
length.  There  was  little  difference  in  the  output  beam  quality;  all  cases 
giving  a  beam  quality  of  <1.02.  Preliminary  calculations  were  also  per¬ 
formed  using  a  beam  propagated  through  the  optical  train  as  input  to  the 
amplifier.  These  calculations  also  showed  a  small  change  in  beam  quality. 

It  is  known  that  large  phase  variations  that  produce  intensity  varia¬ 
tions  in  a  saturable  gain  medium  will  not  be  completely  corrected  by  phase 
conjugation.  Instead  of  the  usual  plane  wave  input  to  LFCM,  a  field  was 
read  in  containing  a  180-degree  phase  jump,  which  was  meant  to  simulate  a 
large  piston  error.  For  high  spatial  resolution,  512  transverse  x  points 
were  used.  It  was  immediately  found  that  the  beam  quality  predicted  by  LFCM 
for  the  output  beam  was  very  poor,  on  the  order  of  5.  The  reason  for  this 
is  that  in  the  real  device  the  output  beam  will  go  back  through  the  source 
of  the  piston  error  and  will  be  corrected,  but  the  LFCM  has  no  mathematical 
description  of  this.  Instead,  a  new  calculation  of  beam  quality  needs  to 
be  done  when  piston  errors  are  present.  The  beam  quality  routine  STREHL 
calculates  the  variance  of  the  phase  of  the  output  beam  amplitude  UOUT(J). 
When  piston  errors  are  present,  the  same  calculation  is  performed  with  the 
new  vector  J  UINJ (J) *U0UT (J)  which  has  the  property  that  if  UOUT  is  the 


conjugate  of  UINJ,  the  phase  di scontinuity  is  cancelled.  The  change  works 
well  for  the  simple  input  piston  errors  and  can  be  used  in  the  case  of  more 
realistic  input  beams. 

Another  problem  with  the  Strehl  calculation  of  beam  quality  was  dis¬ 
covered  when  the  new  code  option  of  a  variable  small  signal  gain  was  exer¬ 
cised.  This  option  caused  some  low  intensities  in  the  output  beam  within 
the  geometric  aperture.  A  peculiarity  of  the  phase  angle  calculations  is 
that  low  intensities  appear  to  yield  large  phase  angles.  These  large 
angles  lead  to  a  large  variance  in  phase  in  the  Strehl  calculation  and, 
hence,  poor  beam  quality.  However,  this  result  is  unphysical,  as  regions 
of  low  beam  intensity  should  obviously  carry  less  weight  than  regions  of 
high  intensity.  Therefore,  it  was  decided  to  intensity-weight  the  sums 
used  to  calculate  the  average  phase  and  variance  in  the  Strehl  routine. 

This  simple  process  brought  about  a  dramatic  improvement  in  the  beam 
quality  when  the  small  signal  gain  varies  in  space. 

5.2.6  Phase  Screens  for  the  LFCM  Code 

In  order  to  account  for  phase  changes  in  the  field  propagating  through 
a  medium  with  an  inhomogeneous  index  of  refraction,  the  original  LFCM  code 
was  modified  to  read  phase  screens,  as  input  data,  which  may  be  applied  at 
various  positions  in  the  direction  of  propagation.  A  separate  screen 
generation  code  was  written  to  accept  either  a  numerical  or  functional 
description  of  the  spatial  variation  of  index  of  refraction.  At  each 
position  at  which  a  phase  screen  is  required,  a  set  of  points  identical  to 
the  sample  points  of  the  LFCM  code  is  specified.  Through  each  po^nt  the 
direction  of  the  ray  from  the  virtual  source  of  the  converging  or  diverging 
field  is  determined  and  the  line  integral  of  (index  of  refraction  -  1.0) 
from  the  given  point  to  the  plane  of  the  next  phase  screen  is  evaluated. 

This  procedure  is  repeated  for  each  sample  point  and  each  phase  screen. 
Special  consideration  was  required  for  ray  paths  which  lie  partially  within 
the  medium  and  partially  outside.  The  calculated  OPD  values  are  then 
written  in  a  format  compatible  with  the  LFCM  code. 


Preliminary  data  values  for  index  of  refraction  have  been  provided  by 
P.  Lohn  (Ref.  5.4  and  5.5),  which  were  the  result  of  previous  flow  calcula¬ 
tions  for  the  ALPHA  configuration.  The  calculated  index  was  averaged  along 
the  z  direction  (parallel  to  the  axis  of  the  cylindrical  amplifier), 
leaving  only  a  dependence  on  radial  and  azimuthal  coordinate.  Since  the 
amplifier  code  models  only  radial  variations,  screens  will  be  constructed 
for  the  two  values  of  azimuthal  coordinate  which  give  the  maximum  and 
minimum  index  values.  The  first  set  of  screens  is  being  constructed. 

5.2.6. 1  Cavity  Medium  Aberrations 

The  preliminary  PALS  design  concept  for  the  amplifiers  embodies  the 
nozzle  designs  and  configuration  used  for  the  ALPHA  program.  Flow-field 
calculations  performed  for  ALPHA  may  be  used  to  provide  an  estimate  of  the 
cavity  medium  index  of  refraction  inhomogeneities.  The  results  are  then  to 
be  used  to  generate  phase  screens  for  the  LFCM  amplifier  code  in  order  to 
model  the  effect  of  the  inhomogeneities  on  field  propagation  through  the 
amplifier  and  on  conjugation  fidelity. 

A  report  on  the  ALPHA  calculations  is  given  in  reference  5.4.  The 
results  are  summarized  below. 

The  calculations  are  expressed  as  OPD  variations  attributable  to  each 
of  numerous  causes  of  density  variations.  If  the  l'th  aberration  produces 
an  index  variation  of  Anj  over  a  distance  ALi  in  the  z  direction  (the  z 
axis  being  parallel  to  the  axis  of  the  cylindrical  amplifier),  then  the  OPD 
due  to  that  source  is,  for  convenience,  written  in  the  form 


0PDi  =  L  (n-1)  A. 


where 


>.4  "ALPHA  OPD  Map  for  PALS,"  P.D.  Lohn,  TRW  IOC  86. K514. 5-009,  12  Feb  1986 
>.5  ALPHA  I  Cavity  Medium  Aberrations,"  P.D.  Lohn,  TRW  IOC  81.4314.3-046 
28  Jul  1981 .  ’ 


and  n  is  the  index  of  refraction  averaged  over  the  length  L  of  the  device. 
Values  of  n-1  are  shown  in  Figure  5-12  as  a  function  of  r,  which  is  the 
radial  coordinate  measured  outward  from  the  nozzle  exit  plane.  Values  of 
Ai  for  a  number  of  aberration  sources  are  listed  in  Table  5-2.  The  nozzles 
are  arranged  in  rings,  with  approximately  13  rings  per  meter  of  amplifier 
length. 

The  number  of  rings,  which  appear  in  the  last  column  of  several 
entries  in  Table  5-2,  are  therefore  determined  by  the  length  of  amplifier 
specified.  The  angle  ^  is  the  polar  coordinate  in  cylindrical  geometry. 
From  Figure  5-12  and  Table  5-2  and  assuming  a  radial  beam  width  of  2-3  cm 
and  a  gain  length  of  9  meters,  an  OPD  difference  across  the  beam  of  ~10~3  A 
is  estimated  for  a  complete  trip  through  the  amplifier,  with  a  maximum  OPD 
of  about  twice  this  amount  for  rays  at  different  azimuthal  positions.  A 
potentially  much  larger  OPD  variation  may  arise  from  the  ALPHA  nozzle  base 
purge  blockage,  discussed  in  reference  5.6.  This  contribution  to  the 
amplifier  medium  aberrations,  and  the  manner  in  which  it  may  be  incorpo¬ 
rated  into  the  amplifier  model,  will  be  examined  in  detail  in  the  next 
phase  of  the  program. 


Figure  5-12.  Average  cavity  index  of  refraction. 


5.6  "Base  Purge  Blockage  OPD,"  G.M.  Harpole,  TRW  IOC  85.K514.3-04] 
28  Jun  1985. 


TABLE  5-2.  GGA  Cavity  Medium  Aberrations. 


Source 

Disturbance 

Symmetry 

or 

Correlation 

Frequency 

(Hz) 

Magnitude 
(An/n-1) AL/L* 

Primary  nozzle 
throat  size 
variations  due 
to  manufactur¬ 
ing  tolerances 

Azimuthal 
mass  flow 
non¬ 
uni  formi ti es 

Uncorrel ated 

0 

0.013  cos4^ 

>J  brings 

Primary  nozzle 
throat  size 
and  composition 
variations: 
thermal ly 
induced 

Azimuthal 
mass  flow 
non- 

uni formi ti  es 

Correlated 

0 

1/2  (0.00136) 
cos20  and 

1/2  (0.00327) 
cos4^ 

Combustor 
noise  based 
on  (VAMP 
results) 

Turbulence 

4-400 

<0.005  (rms) 

Combustor 
acoustic 
osci 1 lation 
based  on 
(VAMP  results) 

Natural 

acoustic 

modes 

<0.005  (rms) 

Secondary 
blade 
injection 
angl  e 
tolerance 

Period  1/6-inch 
uncorrelated 

0 

0.01  (rms) 

y#ri ngs 

Secondary 
blade  blade- 
to-blade  fuel 
flow  variation 

Period  1/6-inch 
uncorrelated 

0 

0.01  (rms) 

v4rings 

Secondary 
blade  flow 
vari ation 
within  blade 

Period  1/6-inch 
uncorrelated 

0 

0.007  (rms) 

\Arings 

Cavi ty 

isolation 

system 

Shear  layers 

4  x  lO"3 

P-V 

*1-2  meters. 


5. 2. 6. 2  Phase  Screen  Calculations  for  Propagation  in  Amplifier 


Electromagnetic  wave  propagation  in  weakly  diffracting  media  may  be 
approximated  by  employing  the  concept  of  geometric  phase  screens.  This 
section  describes  the  computer  code  that  was  developed  to  implement  this 
concept.  Four  or  more  screens  are  typically  used  in  the  LFCM  context  for 
amplifier  calculations. 

5.2.6. 3  Description  of  Phase  Screen  Geometry  Used  for  Calculations 

The  wave(s)  are  assumed  to  emanate  from  an  apparent  source  located  at 
(XS,ZS).  The  medium  is  contained  within  the  boundaries  delimited  by  (XMIN, 
ZMIN) ,  (XMIN,  ZMAX) ,  (XMAX,  ZMIN) ,  and  (XMAX,  ZMAX)  (Figure  5-13 ) . 

The  points  labeled  by  a  square  (□)  are  locations  whose  coordinates  and 
values  of  refractive  index  are  known.  Only  four  phase  screens  are  used  in 
this  example.  The  phase  screens  are  to  be  located  at  the  values  of  Z  coor¬ 
dinates  specified  by  SCRNZS(l),  SCRNZS(2),  SCRNZS(3),  and  SCRNZS(4).  Note 
that  a  phase  screen  may  be  located  outside  the  medium.  The  values  of  the 
x-coordinate  for  which  the  phase  coefficients  are  to  be  computed  are  shown 
as  dark  circles  (#) .  The  line  PO  joining  the  source  point  (XS,  ZS)  with  a 
point  on  the  phase  screen  is  the  ray  for  which  the  phase  coefficient  is 
sought.  The  values  of  the  coefficients  obtained  can  then  be  used  to  propa¬ 
gate  the  field  to  the  next  phase  screen  [SCRNZS  (2)  in  this  case]  and  the 
procedure  repeated  to  calculate  the  phase  coefficients  for  propagation  from 
this  location  to  the  next  phase  screen  at  the  next  position.  The  set  of 
values  of  the  x-coordinate  is  the  same  for  all  screen  positions.  Note  that 
the  rays  may  come  from  apparent  sources  at  various  locations,  such  as  rays 
QA  and  RB.  The  location  of  the  source  (XS,  ZS)  may  be  anywhere  outside  the 
medium. 

The  propagation  of  waves  from  one  phase  screen  position  to  a  phase 
screen  at  another  position  may  engender  several  cases.  These  are  illus¬ 
trated  in  Figure  5-14a-c.  In  each  case  the  medium  is  denoted  in  simplified 
form  by  the  rectangle  and  S  indicates  the  screen  phase  position  in  the 
system. 
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Figure  5-14.  Three  principal  cases  of  phase  factor  calculations  for 
propagation  from  a  point  outside  the  medium  into  the 


medium. 


3 


For  purposes  of  discussion,  only  apparent  sources  from  the  left  are 
considered.  Screen  S]  in  Figure  5-14a  and  screen  S2  in  Figure  5- 14b ,  are 
then,  respectively,  the  first  phase  screens  encountered  for  each  configura¬ 
tion.  In  the  case  of  Figure  5-14a,  the  phase  coefficient  is  to  be  calcu¬ 
lated  for  a  phase  screen  point  outside  the  medium  but  the  wave  propagates 
to  a  screen  point  within  the  medium.  The  partially  broken  line,  continued 
by  the  solid  line,  indicates  the  ray  direction.  Only  the  solid  portion  is 
used  for  physical  length  calculation.  The  refractive  index  is  taken  as  the 
value  at  the  point  on  the  solid  line  marked  by  an  X. 

In  Figure  5-14b,  the  phase  coefficient  is  to  be  calculated  for  a 
screen  located  within  the  medium  and  the  ray  also  propagates  to  a  phase 
screen  within  the  medium.  The  physical  distance  is  given  by  the  solid  line 
and  the  refractive  index  taken  again  to  be  the  value  at  the  point  marked  by 
an  X. 

Figure  5-14c  shows  the  wave  propagating  to  a  screen  outside  the 
medium.  As  before,  X  marks  the  point  for  which  the  refractive  index  is 
used  and  the  physical  path  is  the  solid  line. 

Figure  5-15  is  a  simplified  sketch  of  the  phase  screens  and  a  given 
ray  to  be  used  for  describing  the  computational  procedure.  Three  screens 
are  used  for  illustrative  purposes.  The  direction  of  the  ray  is  given  by 
the  slope  of  the  line  joining  (XS,  ZS)  to  (a,  Sj).  The  latter  phase  screen 
point  is  the  one  from  which  the  field  propagates.  The  point  (Xlf  S2)  is 
the  intersection  of  the  extended  ray  with  the  screen  S2.  This  is  the  phase 
screen  point  to  which  the  wave  propagates.  For  this  particular  example, 
the  phase  factor  is  computed  as 

0] (a)  =  [n ( x )  -  1] .AS 

where  0j(a)  is  the  value  of  the  phase  coefficient  on  screen  Si  at  the  loca¬ 
tion  (a)  and  the  value  of  refractive  index  n(x)  is  the  value  at  the  point 


Figure  5-15.  Sketch  illustrating  computational  procedure  for 
phase  factor  (^(a)  =  [n(x)-l]As.) . 

marked  by  an  x.  In  this  case  AS  is  the  physical  distance  between  the  two 
screens  used  in  the  phase  factor  calculation. 

When  the  inhomogeneity  of  the  medium  varies  rapidly  it  may  be  neces¬ 
sary  to  locate  more  phase  screens  in  the  propagation  path.  Consider  Figure 
5-16;  the  three  previous  screens  are  illustrated  except  now  the  source  is 
assumed  to  be  at  a  different  location.  Because  more  accuracy  is  required, 
additional  intermediate  phase  screens  have  to  be  used  in  the  calculation. 

In  the  figure,  two  intermediate  screens  are  used.  These  are  labeled  aa  and 
ap.  The  phase  coefficient  for  this  example  will  be  computed  as 

N-l 

Ma)  =  H  (n(p)  -  1)  ASp  (5.2) 

P=1  X 


where  ^i(a)  is  the  value  of  the  phase  coefficient  on  phase  screen  S]  at  the 
point  a;  n(p)  is  the  value  of  the  refractive  index  at  the  pth  intersection 
of  the  phase  screen  and  the  ray.  ASpx  is  the  physical  distance  between  the 


point  on  the  pth  screen  to  the  point  on  the  (p+l)st  screen.  N  is  the  total 
number  of  screens;  in  this  case,  N  =  4.  The  additional  interpolating 
screens  are  evenly  spaced  between  the  given  screens  for  which  the  phase 
coefficients  are  to  be  computed. 

5.2.7  Extraction  Efficiency 

This  section  summarizes  the  results  of  an  APACHE-oriented  extraction 
efficiency  study.  These  results  were  obtained  using  a  line-center  (Voigt 
profile)  model.  The  codes  used  to  obtain  the  results,  ANOM  and  CLAM,  and 
their  anchoring,  are  also  described. 

5.2. 7.1  ANOM  Code 

The  ANOM  code,  developed  under  the  Advanced  Chemical  Laser  Optics 
Study  (ACLOS)  program,  is  a  one-dimensional  oscillator  model.  It  follows 
intensity  and  predicts  gain  along  a  one-dimensional  path  through  the 
oscillator,  modeling  the  feedback  process  with  a  reflectivity  factor.  Two 
cascade-connected  transitions  are  modeled.  The  gain  model  predicts  v  =  2, 

1,  and  0  vibrational  level  populations  (N2,  Nj ,  and  Nq)  via  a  kinetics 
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in  which  the  I's  are  intensities,  Aug  is  the  Doppler-broadened  linewidth, 
B ( V , J )  is  the  Einstein  B  coefficient  for  transition  Py+i(J),  and 


Bj  =  (he  Be/kT)  exp  [-J(J+1)  he  Be/kT] 


the  equilibrium  statistical  weight  of  rotational  state  with  quantum  number  J 

The  ANOM  gain  model  was  anchored  to  CROQ  code  results  at  a  flow 
coordinate  of  1.3  cm.  The  F,  R,  and  Q  parameters  were  adjusted  until  the 
ANOM  code  vibrational  level  populations  Nq,  Nj,  and  N2  matched  those  pre¬ 
dicted  by  the  CROQ  code.  (CROQ  is  a  three-dimensional  annular  laser  model, 
whose  gain  model  has  a  detailed  kinetics  and  a  simple  flow  model.)  Fitting 
was  done  at  gain  lengths  of  200,  518,  and  916  cm.  The  parameters  are  a 
function  of  gain  length  (Gl,  cm)  as  follows: 

F0  =  °* 

F]  -  1.0  x  1020 

F2  =  7.0  x  1020 

R0  =  (8.19  -  4.73x10-3  GL)  x  10^0 

Rj  =  (4.81  +  4.73xlO-3  GL)  x  1020 

R2  =  11.0  x  1020 

Q0  =  3.0  x  105 

Ql  =  (2.94  +  7.77x10-3  GL)  x  10$ 

QZ  =  (13.17  -  1.66x10-16  GL  -  1 . 74xl0"5  GL2)  x  105. 

Further  details  of  the  gain  model  (e.g.,  how  rotational  and  velocity 
level  populations  and  gain  are  calculated)  are  given  in  the  ACL0S  final 
report  (Ref.  5.7) 


5.7  "Advanced  Chemical  Laser  Optics  Study  (ACL0S),"  Final  Report, 
AFWL-TR-82-54 . 


Oscillator  Model 


ANOM  models  a  one-dimensional  path  through  the  oscillator.  There  is 
one  “gain  sheet."  For  a  standing-wave  resonator  (as  opposed  to  a  ring 
resonator),  the  intensity  at  the  back  mirror  is  multiplied  by  exp(G  CL)  (CL 
being  the  cavity  length  and  G  the  gain)  to  "propagate"  to  the  scraper;  this 
intensity  is  then  multiplied  by  R,  the  feedback  ratio,  and  again  multiplied 
by  exp(G  CL)  to  "propagate"  back  to  the  back  mirror.  The  resulting  inten¬ 
sity  is  compared  to  the  starting  intensity  for  self-consistency;  if  they 
are  not  equal  to  within  the  convergence  criterion,  the  intensity  is  updated 
according  to: 


I 


new 


(1-d) 


[(!old)  exp(G  CL)  R  exp(G  CL)]d 


where  d  is  a  user-specified  damping  factor.  The  intensity  used  to  calcu¬ 
late  the  gain  G  is  the  sum  of  the  average  intensities  for  each  pass: 


1avg(pass  1}  =  1  [exp(G  CL)  -  1]  /  (G  CL) 
and 

Iavg(pass  2)  =  I  exp (G  CL)  R  [exp(G-CL)  -  1]  /  (G  CL). 

Formulae  for  the  gain  calculation  are  given  in  reference  5.7.  When  the 
intensity  at  the  back  mirror  has  converged,  the  output  intensity  at  the 
scraper  is  calculated  as: 


lout  r  1  e>p(G  CL)  (1  -  R). 

Although  the  ALPHA  resonator  is  a  ring,  it  is  modeled  as  a  standing 
wave  resonator  because  the  photons  pass  through  the  gain  medium  in  both 
directions.  The  resonator  length  is  input  to  the  ANOM  code  as  half  of  its 
actual  value,  so  that  the  proper  longitudinal  mode  spacing  for  the  ring  is 
simulated.  Longitudinal  mode  spacing  is  approximated  as  that  of  a  "closed 
cavity"  resonator. 

The  ANOM  code  output  intensities  are  used  as  input  to  the  amplifier 
model,  CLAM.  A  mul t ipl e- longi tud i nal -mode  output  from  an  ANOM  run 
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simulating  a  200-cm  gain  length  ALPHA  resonator  is  given  in  Figures  5-17a 
through  5-17d. 

5. 2. 7. 2  CLAM  Code 


CLAM  predicts  the  amplification  of  the  oscillator  intensity  by  the 
amplifier.  The  CLAM  gain  model  is  the  same  as  that  of  the  ANOM  code. 


Amplifier  Model 


The  CLAM  amplifier  model  is  an  adaptation  of  the  amplifier  model  in 
the  code  PSTAR.  It  was  modified  so  that  it  would  model  an  ALPHA-1  ike 
device.  A  schematic  of  this  device  is  shown  in  Figure  5-18. 


If  the  logical  parameter  FORWARD  is  true,  the  beam  from  the  oscillator 
enters  at  the  waxicon  inner  cone,  is  annularized  to  the  waxicon  outer  cone, 
then  expands  as  it  passes  through  the  gain  medium  to  the  rear  element;  it 
then  gets  smaller  as  it  passes  back  through  the  gain  medium  to  the 
reflaxicon  outer  cone,  is  compacted  to  the  reflaxicon  inner  cone,  and  is 
sent  to  the  phase  conjugation  device;  upon  return  from  the  phase  conjugator 
(PC),  the  beam  retraces  itself  and  exits  the  amplifier  from  the  waxicon 
inner  cone.  If  FORWARD  is  false,  the  beam  enters  and  exits  at  the 
reflaxicon  inner  cone,  and  the  phase  conjugator  is  on  the  "waxicon  side." 
Note  that  an  actual  ray  in  this  amplifier  would  have  its  upstream  and  down¬ 
stream  passes  through  the  gain  medium  on  opposite  sides  of  the  gain  gener¬ 
ator.  However,  since  our  model  follows  a  one-dimensional  path  through  the 
amplifier,  and  since  it  is  assumed  that  the  intensity  is  cyl i ndrical ly 
symmetric,  the  intensity  path  in  CLAM  will  be  the  one  shown  in  the  figure. 


The  logical  variable  ALPHA2  (for  2  passes  through  the  gain  medium  from 
input  to  the  PC)  controls  whether  or  not  this  model  is  used.  There  is  also 
an  ALPHA1  (1  pass)  option,  which  would  model  the  case  of  a  device  with  one 
axicon  on  either  end  of  the  annular  leg.  This  option  was  rejected  for  the 
ALPHA  resonator  device  because  of  the  extreme  sensitivity  to  axicon  tip-to- 
tip  alignment;  when  the  axicon  inner  cones  are  back-to-back  and  the  outer 
cones  are  in  one  piece,  this  sensitivity  is  ameliorated.  It  may  be,  how¬ 
ever,  that  the  amplifier  would  not  have  such  extreme  sensitivity,  and  that 
the  1-pass  conf iguration  may  have  a  greater  extraction  efficiency.  It 
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Figure  5-17.  ANOM  multi  pi  e-longitudinal -mode  output. 


is  recommended  that  1-pass  extraction  efficiency,  alignment  sensitivity, 
and  polarization  implications  be  investigated  by  the  APACHE  program. 

The  input  intensities  are  scaled  from  the  compact  beam  size  to  the 
annular  beam  size  upon  entering  the  annular  region,  and  scaled  conversely 
when  leaving  the  annular  region.  Furthermore,  as  the  beam  passes  through 
the  annular  region,  the  intensities  and  beam  overlap  are  scaled  according 
to  the  area  change  and  the  proper  geometry,  as  shown  in  the  figure.  In  the 
runs  made  to  date,  a  compact  beam  radius  of  12.5  cm  was  used,  in  order  to 
have  a  uniform  set  of  intensity  results.  However,  in  "real  life"  the 
axicon  inner  cone  sizes  will  be  dictated  by  engineering  considerations.  In 
the  ALPHA2  configuration,  the  beam  sizes  have  been:  rear  element,  2.6  cm; 
waxicon  outer  cone,  1.0  cm;  reflaxicon  outer  cone,  1.3  cm  (as  in  the  ALPHA 
resonator).  It  is  assumed  that  these  will  be  the  best  values  for  use  with 
the  nominal  ALPHA  flow  conditions;  if  the  flow  is  changed,  these  values 
will  have  to  be  reoptimized. 
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There  is  an  arbitrary  (user-specified)  number  of  “gain  sheets."  The 
gain  at  each  "sheet"  is  calculated  using  the  average  intensity  there.  This 
average  intensity  is  made  up  of  contributions  from  the  appropriate 
intensity  going  to  and  coming  from  the  PC  (the  "two-way  intensity"),  and 
the  overlap  factor  times  the  "two-way  intensity"  of  the  opposing  annular 
leg  beam.  As  in  ANOM,  the  average  intensity  is 

lave  =  1  [exp(G  DL)  -  1] / (G  DL) 

where  DL  is  the  gain  length  represented  by  one  gain  sheet. 

Because  the  gain  at  a  particular  gain  sheet  is  determined  by  all  of 
the  intensities  passing  through  the  gain  sheet,  there  needs  to  be  an 
initial  guess  for  all  of  the  intensities.  CLAM  subroutine  GUESS  supplies 
these  initial  intensity  guesses.  GUESS  solves  a  one-gain-sheet  amplifier 
(in  a  way  analogous  to  that  in  which  the  ANOM  code  solves  a  one-gain-sheet 
oscillator),  then  defines  the  intensity  at  gain  sheet  i  to  be 

Ij  =  I , exp(G  DL). 

The  PC  is  modeled  by  a  reflectivity,  Rpo  Thus 

I  (after  PC)  =  I  (before  PC)  x  RpQ. 

The  CLAM  code  also  has  a  logical  parameter,  ONEWAY,  which  allows 
modeling  of  a  "one-way"  pass  through  the  amplifier;  i . e . ,  the  beam  enters 
at  one  axicon,  goes  through  the  gain  medium,  then  leaves  at  the  other 
axicon.  This  allows  the  code  to  simulate  one  pass  through  an  oscillator; 
this  oscillator  simulation  is  used  in  the  intensi ty-to-power  calibration 
described  below. 

Intensity/Power  Calibration  and  Consistency  Checks 

a.  ANOM/CLAM  Consistency  Check 

The  ALPHA  device  has  an  outcoupling  fraction  of  0.88  and  a  feedback 
fraction  of  0.12  (as  predicted  by  the  CROQ  code).  Thus  the  ANOM  code  was 
run  with  R  =  0.12  in  order  to  model  the  ALPHA  resonator.  Since  the  ANOM 
code  output  assumes  the  0.88  outcoupling  fraction,  if  the  ANOM  intensity 
output  is  multiplied  by  0.12/0.88  =  0.136,  that  intensity  should  represent 
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the  amount  of  intensity  fed  back  into  the  oscillator  on  each  pass.  If  the 
CLAM  model  is  configured  as  an  oscillator  (ONEWAY  true)  with  the  same  gain 
length,  then  this  input  intensity  should  be  amplified  by  a  factor  of 
1.0/0.12  =  8.33,  the  same  as  the  one-pass  amplification  of  the  oscillator. 

Since  the  ANOM  code  does  not  scale  the  annular  intensities  to  the 
compact  beam  area,  the  CLAM  code  was  setup  to  give  the  same  area  for  both 
the  compact  beam  and  the  beam  entering  the  annular  leg.  The  lines  modeled 
were  P 2(7)  and  Pi(8).  Table  5-3  gives  the  results  of  the  CLAM  runs  for 
both  line  center  and  multiple  longitudinal  mode  models.  The  probable 
reason  that  the  results  are  not  exactly  8.33  is  that  ANOM  does  not  have  the 
ALPHA  device  annular  leg  geometry. 

TABLE  5-3.  ANOM/CLAM  Consistency  Check. 


Line 

Line  Center  Model 

Multiple  Mode  Model 

<W!1n 

should  be  8.33) 

P2(7) 

8.32 

8.38 

Pi  (8) 

8.29 

8.35 

b.  ANOM/CRQQ  Consistency  Check 

Since  the  ANOM  "kinetics  parameters"  have  been  adjusted  so  that  the 
ANOM-predicted  values  of  Nq,  N^,  and  N2  match  those  of  CROQ,  the  consis¬ 
tency  check  must  involve  the  predicted  intensities. 

The  CROQ  code  predicts  that  the  lines  P2 ( 7 )  and  Pj(8)  are  two  of  the 
most  powerful  lines  in  the  resonator  output  for  gain  lengths  ranging  from 
200  to  over  900  cm.  Therefore,  these  lines  were  modeled  with  ANOM  and 
CLAM.  ANOM  calls  the  upper  level  intensity  II  and  the  lower,  10.  If  ANOM 
correctly  models  the  gain  medium  response,  then  11/10  is  expected  to  be  the 
same  as  that  of  CROQ  and  11+10  to  be  proportional  to  that  of  CROQ  (with  the 
same  proport iona 1 i ty  constant  for  all  gain  lengths). 

Table  5-4  gives  the  results  of  these  tests,  and  shows  the  details  of 
the  anchoring  results  for  vibrational  level  populations.  The  intensity 
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TABLE  5-4.  ANOM/CROQ  Consistency  Check. 


Gain  Length  (cm) 

204 

(ALPHA  I) 

518 

(ALPHAII) 

705 

("APACHE7") 

916 

("APACH 

Gain  Model  Anchoring 

N0  (CROQ)  x  10-15 

3.1 

3.3 

3.4 

Nq  (ANOM)  x  10-15 

3.0 

3.2 

3.4 

Ni  (CROQ)  x  10-15 

1.9 

1.8 

1.8 

Nj  (ANOM)  x  10-15 

1.9 

1.8 

1.8 

N2  (CROQ)  x  10-15 

1.1 

1.0 

1.0 

N2  (ANOM)  x  10"15 

1.2 

1.1 

1.0 

Intensity  Ratio 
11/10  (CROQ) 
II / 10  (ANOM) 


PCROO^Ann’  Area^ 
11+10  (ANOM) 


results  are  only  fair.  The  result  at  the  200  cm  gain  length  is  somewhat 
different  from  those  at  the  other  gain  lengths  due  to  the  much  larger  power 
fraction  on  P2(7)  --  33°,  at  200  cm;  22°,  to  241  at  the  longer  gain  lengths. 

c.  Scalinq  ANOM-Predi cted  Intensity  to  CROQ  Power 


The  ANOM  code  models  only  the  annular  leg  of  the  ALt'HA  resonator.  It 
does  not  model  the  actual  annular  leg  geometry,  but  assumes  total  overlap 
of  both  passes  through  the  gain  medium,  and  a  constant  beam  width.  If  the 
annular  leg  were  modeled  properly  in  ANOM,  the  annular  output  intensity  to 
the  compact  intensity  could  be  simply  scaled  by  the  ratio  of  the  areas. 
However,  there  is  no  way  of  knowing  the  "effective  gain  width"  (Weff)  for 
the  ANOM  run. 


I, 


If  Weff  were  known,  then  the  ALPHA  I  device  (gain  length  200  cm)  out¬ 
put  power  P}  would  be  represented  by  the  intensity 

(II  +  10)  (from  ANOM)  x  (AREAannuiar) / (AREAcompact) 

where 

AREAannuiar  =  2  t  Ravg  Weff 

Ravg  =  R  (gain  generator)  +  Mode  Offset  +  1/2  Weff 
=  approximately  56.6  for  ALPHA  I;  thus 

AREAannular  =  356  Weff 

AREAcompact  =  *  Rcompact^  =  *  (12.5)2  =  491, 

Since  Weff  could  reasonably  be  anything  between  1.3  and  2.6,  the  ratio 
AREAannu]ar/AREAconipact  is  expected  to  be  between  0.94  and  1.89. 

In  order  to  actually  scale  the  intensities,  the  fact  that  the  one-pass 
amplification  of  the  oscillator  is  8.33  was  used.  The  oscillator  inten¬ 
sities  were  adjusted  (by  multiplying  by  the  scaling  factor  IMULT)  until  the 
CLAM  output  intensities  for  the  200  cm  gain  length  (ALPHA2  true,  FORWARD 
true,  ONEWAY  true)  were  8.33  times  the  input  intensities.  The  IMULT  value 
needed  to  give  this  result  was  0.15,  very  close  to  the  expected  value  of 
0.136  (see  AN0M/CLAM  Consistency  Check,  above).  The  CLAM  output  intensity 
sum  represents  1.0/0.88  of  Pj ,  since  CLAM  does  not  subtract  any  of  its 
output  for  feedback  as  ANOM  does.  Thus 

(II  +  10)  (CLAM  output)  x  0.88  represents  Pj. 

The  actual  numbers  were: 

(II  +  10)  (CLAM  output)  x  0.88  =  1889 


(II  +  10)  (ANOM  output) 
Ratio  (Scaling  Factor) 


=  1136 


=  1.66 


which  gives  a  reasonable  number  for  Weff,  2.3  cm. 

Thus,  in  order  to  simulate  an  amplifier  input  power  P-j  =  f  Pj,  the 
ALPHAI  ANOM  output  should  be  scaled  by  IMULT  =  f  1.66;  the  output  power  P0 
is  determined  by 
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P0  =  (II  +  10)  (CLAM  output)/1889  x  Pj. 

Extraction  Efficiency  for  Line  Center  (Voigt  Profile)  Model 

a.  Amplifier  Size  Choices 

The  following  assumptions  were  made: 

1.  The  goal  for  power  output  per  amplifier  is  Pa  =  4.1667  P2, 
where  P2  is  the  ALPHA  II  output  power  goal. 

2.  The  "nominal  ALPHA  flow  conditions"  are  used. 

3.  The  power  predicted  by  CROQ  for  an  oscillator  is  approximately 
the  same  as  the  power  output  of  an  amplifier. 

With  these  assumptions,  the  curve  for  the  required  gain  generator  radius  as 
a  function  of  gain  generator  length  is  shown  in  Figure  5-19. 

The  actual  curve  will  be  similar  to  this  one.  If  the  extraction  is 
more  efficient  than  this  curve  predicts,  then  the  actual  curve  will  be 
lower  and  to  the  left  of  the  curve  shown. 

There  are  a  number  of  factors  which  put  an  upper  limit  on  the  gain 
generator  radius.  These  include:  vehicle  size  limitation,  mirror  manu¬ 
facture  capability,  and  alignment  tolerance.  The  vehicle  size  limitation 
on  gain  generator  radius  is  about  150  cm. 

Also,  a  number  of  factors  put  an  upper  limit  on  the  gain  generator 
length.  These  are  amplified  spontaneous  emission  (ASE),  gain  generator 
structural  stability  and  alignment  tolerance  (which  get  worse  as  the 
length-to-diameter  ratio  gets  larger),  diffraction  loss,  the  "gain-phase 
effect,"  and  vehicle  size  limitation.  The  upper  limit  on  gain  generator 
length  is  not  yet  clear. 


Because  the  CROQ  code  has  only  one  gain  sheet,  the  code  is  not  very 
stable  at  gain  lengths  above  9  m.  For  that  and  the  above  reasons, 
amplifier  extraction  efficiency  in  the  gain  generator  length  range  of  up  to 
9  m  was  chosen  for  study. 

Also  it  was  assumed  that  the  annular  leg  mirror  configuration  is  the 
same  as  that  in  the  ALPHA  device,  since  that  configuration  was  shown  to 
produce  the  best  extraction.  Thus  the  geometric  mode  is  assumed  to  be 
2.6  cm  wide  and  be  offset  from  the  gain  generator  by  3  mm.  The  waxicon 


outer  cone  is  assumed  to  be  1.0  cm  wide,  and  the  reflaxicon  outer  cone  is 
assumed  to  be  1.3  cm  wide. 


ALPHA  I-Sized  Amplifier  --  Comparison  with  gn-I<;at  Model 


Before  proceeding  to  the  APACHE-sized  device  studies,  CLAM  was  used  to 
duplicate  the  PALS  MOPA  system  model  (Ref  5.8).  This  study  looked  at  power 
to  the  phase  conjugator  as  a  function  of  power  input  to  the  amplifier. 


The  model  assumed 
g  =  g0/(l  4  I/Isat), 

where  I  is  the  two-way  intensity.  The  model  originally  did  not  account  for 
overlap  of  the  annular  leg  passes,  however  it  was  subsequently  modified  and 
the  new  results  are  presented  in  Section  5.2.8. 

The  CLAM  calculation  was  done  both  with  and  without  annular  leg  over¬ 
lap.  Figure  5-20  shows  the  results.  (The  numbers  by  the  symbols  are  the 


5.8  "PALS  MOPA  System  Model,”  V.W.  Chai,  TRW  IOC  86 . K325 , VWC-001 , 
15  Jan  1986. 
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Figure  5-20.  Extraction  for  ALPHA  I-size  amplifier. 

output  powers.  All  powers  have  been  normalized  to  Pi.)  The  results  obey 
the  approximate  algorithm 

PPC  =  (Pin  Pout  /  RPC)1/2 


where  Rpc  is  the  effective  pulse  conjugator  reflectivity. 

It  is  seen  that  the  annular  leg  overlap  leads  to  about  20  percent 
lower  power  to  the  phase  conjugator  than  if  no  overlap  were  present.  Als 
it  appears  that  g0  and  Isat  values  of  0.025  and  700,  respectively,  would 
give  gain  model  results  similar  to  those  of  CLAM.  (These  are  lower  than 
the  values  of  0.03  and  1100  based  on  numbers  predicted  by  the  CR0Q  code.) 

c.  Extraction  Efficiency  at  5,  7,  and  9  m  Gain  Lengths 


Figure  5-21  shows  the  line-center-model  predictions  of  CLAM  for  the 
power  out  and  power  to  the  phase  conjugator  of  5,  7,  and  9  m  amplifiers. 


HS 


Figure  5-21.  Ppr  and  PmiT  as  a  function  of  Pjn  (CLAM,  Voigt  profile 
model 7^ 

The  gain  generator  area  was  kept  constant.  Powers  shown  are  normalized  to 
P/\  =  4.1667  ?2-  The  9ain  length  does  not  seem  to  have  a  significant  effect 
on  extraction  efficiency. 

The  output  power  appears  to  be  a  weak  function  of  input  power.  Thus, 
the  requirement  for  input  power  will  be  determined  by  the  requirement  for 
power  to  the  phase  conjugator  and  possibly  by  ASE  considerations. 

5.2.8  PALS  MOPA  System  Model  Upgrade 

The  PALS  MOPA  system  model  documented  in  reference  5.8  is  a  simple 
microcomputer  model  that  has  been  modified  to  incorporate  the  beam  overlap 
as  the  beam  travels  back  and  forth  through  the  ALPHA-1  ike  gain  generator. 
This  section  describes  the  changes  made  to  the  model  and  shows  the  effects 
of  the  beam  overlap  to  the  power  delivered  to  the  SBS  cell  and  the  system 
output  power. 

5. 2.8.1  Background 

The  model  documented  in  reference  5.8  assumed  that  there  is  no  overlap 
between  the  annular  leg  passes  through  the  gain  generator,  for  a  first 
order  assessment  of  the  PALS  system  performance.  This  led  to  an  optimistic 
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prediction  of  the  power  produced  by  the  amplifier.  The  effort  described 
here  includes  the  beam  overlap  to  obtain  a  better  estimation  of  the  power 
out  of  the  amplifier.  The  upgraded  model  also  simulates  the  beam  size 
(which  was  assumed  to  be  constant  in  the  previous  effort),  as  it  propagates 
from  the  waxicon  outer  cone  to  the  rear  cone  and  back  to  the  reflexicon 
outer  cone,  using  the  ALPHA  configuration. 

5. 2. 8. 2  Simulation  of  the  Beam  Size 

The  size  of  the  annular  leg  determines  the  effective  annular  area  and 
thus  the  intensity  in  the  gain  medium.  Figure  5-22  shows  the  mirror 
configuration  used  in  the  model  with  a  variable  gain  generator  length.  The 
beam  expands  linearly  from  1.0  cm  at  the  waxicon  outer  cone  to  2.6  cm  at 
the  rear  cone  mirror  and  contracts  linearly  from  2.6  to  1.3  cm  at  the 
reflaxicon  outer  cone.  The  model  ignores  the  curvatures  of  the  mirrors 
since  they  have  a  negligible  effect  on  the  total  length.  A  typical  example 
of  the  annular  area  is  shown  in  Figure  5-23  as  a  function  of  the  position 
in  the  gain  generator. 

5.2.8. 3  Modelinq  of  the  Beam  Overlap 

The  model  first  determines  where  the  overlap  occurs  within  the  gain 
generator.  The  inception  of  the  overlap  is  calculated  when  the  summation  of 
the  annular  thicknesses  of  the  passes  amounts  to  an  arbitrary  value  of  2.7 
cm.  Note  that  the  annular  leg  passes  are  bounded  between  2.6  cm  at  the 
rear  cone  mirror  and  2.8  cm  at  the  waxicon/ref 1 axicon  outer  cone.  The 
overlap  factor  then  varies  linearly  from  zero  at  the  inception  point  to  1.0 
at  the  rear  cone  mirror.  The  calculated  overlap  factor  is  shown  in  Figure 
5-24  as  a  function  of  the  position  in  a  7-meter  gain  generator.  It  is 
shown  that  there  is  no  overlap  at  the  waxicon, reflaxicon  end  of  the  gain 
enerator. 

As  described  in  reference  5.8,  the  equivalent  gain  at  any  position,  x, 
is  calculated  by 

<  -  90'd  *  noT.'isat) 
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Figure  5-24.  Typical  beam  overlap  factor. 

where  x  and  n  denote  the  position  being  calculated  and  the  current  iter¬ 
ation,  respectively.  The  ITOT ,  the  total  intensity  present  at  any  point 
is  calculated  according  to  the  following  algorithms: 

For  the  first  pass,  unconjugated  beam, 


ITOT  =  If  " 
x-1 


X  F 


x 


For  the  second  pass,  unconjugated  beam, 


ITOT  =  If  " 
x-1 


X  F 


x 


For  the  second  pass,  conjugate  beam, 


For  the  first  pass,  conjugate  beam, 


where  F  represents  the  overlap  factor  and  the  subscripts  f,  b,  and  p  denote 
the  forward  nonconjugated  beam,  backward  conjugated  beam,  and  beam  pass 
(first  or  second)  respectively. 

The  instantaneous  intensity  at  any  point  is  then  calculated  as 
follows : 

For  the  forward-going,  unconjugated  beams, 

Ifn  =  If  (1  +  g"  Ax)  (areaxl/areax) 

For  the  backward-going,  conjugate  beams, 

]b"  =  (1  +  gx  Ax)  (areax+i/area>> 

5. 2.8.4  Resul ts 

Figure  5-25  shows  the  equivalent  gain  through  the  gain  generator  for 
the  conditions  stated.  The  equivalent  gain  for  the  case  without  beam 
overlapping  was  also  included  for  comparison.  It  is  shown  that  the  gain  is 
lower  with  beam  overlapping  during  the  second  pass  compared  to  that  without 
beam  overlap.  The  lower  gain  translates  to  lower  amplification,  resulting 
in  lower  power  produced  by  the  amplifier. 

Figure  5-26  shows  the  power  delivered  to  the  SBS  cell  as  a  function  of 
the  gain  generator  length  (using  the  space-based  laser  normalization  fac¬ 
tor).  The  beam  overlap  resulted  in  15-22  percent  lower  power  delivered  to 
the  SBS  cell.  This  finding  is  in  good  agreement  with  that  documented  in 


reference  5.9.  Figure  5-27  illustrates  the  system  output  power  as  a  func¬ 
tion  of  the  gain  generator  length.  The  beam  overlap  leads  to  25-35  percent 
lower  system  output  power. 

5. 2. 8. 5  Conclusions 

The  MOPA  system  model  has  been  modified  to  incorporate  the  capability 
to  simulate  the  beam  size  as  it  travels  through  the  gain  generator  and  to 
study  the  effect  of  the  beam  overlap  on  the  power  produced  by  the 
amplifier.  For  the  range  of  gain  generator  length  studied,  the  MOPA  system 
model  predicts  a  15-22  percent  lower  power  delivered  to  the  SBS  cell  and  a 
25-35  percent  lower  system  output  power  with  the  beam  overlap  compared  to 
those  without  beam  overlap. 

5.3  LOW/HIGH  POWER  BEAM  SEPARATION  GRATING  RHOMB  ISOLATION  STUDY 

An  analytical  and  subscale  fabrication  study  was  undertaken  to  assess 
performance  and  fabrication  issues  associated  with  using  a  grating  rhomb 
which,  when  properly  configured,  serves  to  couple  the  low  power  beam  from 
the  APACHE  oscillator  to  the  amplifier/phase  conjugation  subsystem,  extract 
the  return  phase  conjugated  high  power  APACHE  beam,  and  also  isolate  the 
APACHE  oscillator  from  the  high  power  return  beam. 

A  representative  grating  rhomb  configuration  which  illustrates  the 
objective  is  shown  in  Figure  5-28.  A  low-power  oscillator  beam  is  spec¬ 
ularly  reflected  onto  the  secondary  of  the  grating  rhomb  which  possesses 
very  high  efficiency  in  the  TE  first  order  polarization.  The  TE  first 
order  is  then  reflected  off  the  grating  primary  and  propagates  to  the  sys¬ 
tem  amplifier.  The  high-power  phase  conjugated  return  beam  with  its  polar¬ 
ization  rotated  by  90  degrees  is  outcoupled  via  the  high  efficiency  TM 
zeroth-order  polarization  of  the  primary.  The  TM  first  order  is  designed 
to  have  very  low  efficiency  so  that  only  a  small  percentage  of  the  high 
power  beam  is  fed  back  to  the  oscillator  subsystem.  Further  isolation  of 
this  fed-back  beam  with  a  separate  element  may  be  necessary. 


5.9  "Extraction  Efficiency  I,"  D.  Dee,  TRW  IOC  85.K323-78, 
31  Dec.  1985. 


The  primary  analytical  tool  used  in  this  study  was  the  in-house  French 
Grating  Code  which  accepts  input  design  parameters  (wavelength,  incident 
angle,  grating  reflection  coefficient,  groove  depth,  period,  duty  cycle, 
and  contour)  and  predicts  grating  efficiencies  for  possible  orders  in  both 
polarizations. 

The  grating  subscale  effort  was  to  explore  fabrication  methods  by 
which  deep  grooves  could  be  milled  into  subscale  samples  using  the  existing 
ALPHA  grating  master.  Controlling  groove  depth  and  groove  contour  and 
measuring  the  effect  on  TE,  TM,  and  zeroth  order  efficiencies  were  initial 
goals. 


5.3.1  Grating  Rhomb  Design  Goals 

The  design  goals  as  set  forth  in  the  PALS  Preliminary  Reguirement 
Review  (PRR)  of  23  October  1985  were: 

1.  Grating  rhomb  efficiency  for  the  incoupled  oscillator  beam  to  be 
greater  than  90  percent 

2.  Grating  rhomb  efficiency  for  the  outcoupled  high  power  beam  to  be 
greater  than  95  percent 

3.  Grating  rhomb  efficiency  for  the  high  power  orthogonal 
polarization  return  beam  to  be  less  than  5  percent. 

These  design  goals  were  arbitrarily  set  pending  more  mature  system  studies. 

5.3.2  French  Grating  Code  Results 

Early  modeling  of  a  lamellar  (sguare-groove)  grating  with  deep  grooves 
showed : 

1.  At  an  incident  angle  of  40  deg  and  using  a  lamellar  (sguare 

groove)  groove  depth  of  1.32  microns  with  a  groove  period  of  3.2 
microns,  it  was  possible  to  in-couple  81  percent  of  the  oscillator 
beam  power  (Figure  5-29,  TE-first  order)  and  outcouple  81  percent 
of  the  high  power  beam  (four  power  weighted  ALPHA  strong  lines) 
while  allowing  only  0.04  percent  of  the  orthogonal  polarized  high 
power  beam  to  return.  (Figure  5-29  TM-first  order). 

2.  Higher  angles  of  incidence  significantly  reduced  oscillator  in¬ 
coupling  efficiencies  (Figure  5-29). 
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Lamellar  groove  first  order  grating  efficiencies. 
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From  these  results  it  was  necessary  to  address  the  problem  of 
increasing  the  high  power  beam  outcoupled  efficiency.  Study  of  the  code 
outputs  suggested  that  outcoupled  efficiency  could  be  enhanced  by 
decreasing  the  angle  of  incidence.  Also,  earlier  work  for  lamellar 
gratings  plotting  zeroth  order  TM  versus  the  quantity  lambda/groove  period 
suggested  that  optimal  output  coupling  could  be  achieved  for  a  spread  of 
wavelengths  by  keeping  the  quantity  lambda/groove  period  between  0.64  and 
0.71  (the  results  presented  in  the  PALS  PRR  used  lambda/groove  period 
values  between  0.86  and  0.91). 

Therefore,  a  groove  period  value  of  4.1  microns  was  chosen  which 
yields  a  lambda/groove  period  between  0.67  to  0.71  for  the  four  strong 
ALPHA  lines  and  code  runs  for  20-degree  incidence  were  performed.  Figure 
5-30  shows  that  a  groove  depth  of  1.34  microns  yields  a  high-power  out- 
coupling  efficiency  off  the  TM  zeroth  order  of  92  percent,  an  oscillator 
incoupling  efficiency  off  the  TE  first  order  of  79  percent,  and  a  TM  first 
order  return  efficiency  of  0.04  percent. 

It  is  to  be  noted  that  the  efficiency  curves  vary  somewhat  drastically 
with  groove  depth.  This  imposes  strict  fabrication  tolerances  on  unifor¬ 
mity  of  groove  depth  over  the  active  area  of  the  grating. 

In  the  course  of  choosing  a  grating  for  the  CLPCT  Oscillator  Sub¬ 
system,  a  sinusoidal  groove  holographic  grating  (American  Holographic)  was 
found  which  was  purported  to  have  very  high  TE  zeroth  order  efficiency 
(95%)  and  low  TE  first  order  efficiency  on  the  grating  code.  Using  a 
groove  period  of  1.67  microns  (6001/mm)  and  incident  angles  of  55,  60,  and 
62  degrees  (reasonable  grating  rhomb  angles),  the  results  shown  in  Figures 
5-31  through  5-33  were  obtained.  For  clarity,  only  the  Pj(7)  and  P2W 
line  efficiencies  are  plotted;  the  Pj (8)  and  P2 ( 7 )  plots  lie  within  the 
efficiency  bands  shown.  First,  it  is  noticed  that  the  curves  do  not  drop 
or  rise  precipitously  as  a  function  of  groove  depth  as  in  the  case  of  the 
lamellar  curves.  This  is  very  desirable  from  the  standpoint  of  fabrication 
tolerancing.  Second,  the  angular  alignment  tolerance  is  large,  on  the 
order  of  several  degrees.  Third,  shallow  groove  depths  in  these  runs 
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Figure  5-32.  Sinusoidal  groove  first  and  zero  order  efficiencies . 
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enhance  TE  zeroth  order  outcoupled  efficiency  and  minimize  TE  first  order 
efficiency  for  the  return  beam,  but  at  the  expense  of  lower  efficiency  TM 
first  order  efficiency  for  the  incoupled  oscillator  power  (the  remainder  of 
the  oscillator  power  being  diverted  to  the  TM  zeroth  order). 

Table  5-5  shows  expected  rhomb  performance  for  both  lamellar  and 
sinusoidal  groove  gratings.  Values  of  grating  efficiency  have  been  squared 
to  account  for  the  rhomb  configuration.  The  only  exception  is  the  high- 
power  outcoupled  beam  which  outcouples  via  the  zeroth  order  of  primary 
only. 

An  inexpensive  sinusoidal  test  grating  will  be  procured  and  character¬ 
ized  to  verify  the  grating  code  predictions. 

5.3.3  Grating  Subscale  Results 

The  basic  materials  procedure  stated  in  the  PALS  PRR  was  used  to  study 
the  effects  of  argon  ion  milling  bare  molybdenum,  and  nickel  and  gold 
coated  substrates  using  the  ALPHA  master.  The  groove  period  of  the  master 
is  3.2  microns. 

Deep  ion  etching  bare  molybdenum  of  nominal  and  good  stock  revealed  an 
undesirable  surface  crazing.  Previous  work  by  others  in  deep-groove  milling 
suggested  that  there  would  be  processing  problems  with  whisker  formation, 
and  this  was  found  to  be  the  case  in  working  with  ni ckel -coated  substrates. 
Suggestions  on  how  to  remove  the  whiskers  included  ul trasoneration,  ion 
milling  the  sample  at  grazing  angle,  and  stripping  the  surface  with  a  thin 
film  setting  agent  (ex.  collodian).  Finally,  several  efforts  were  made  to 
ion  mill  pure  and  1  percent  cobalt  doped  gold  coated  samples.  In  both 
cases,  it  was  found  that  as  groove  depths  of  about  1  micron  wire  approached, 
a  surface  mottling  appeared  indicating  nonuniform  surface  milling. 

5.3.4  Summary 

The  results  of  the  study  to  assess  performance  and  fabrication  issues 
associated  with  use  of  a  grating  rhomb  are  summarized  below. 

1.  The  grating  efficiency  curves  for  the  lamellar  groove  design  vary 
precipitously  with  groove  depth.  This  sets  tight  fabrication 
tolerances  in  groove  depth  uniformity.  In  addition,  the  subscale 
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TABLE  5-5.  Requirements  Versus  Capabilities. 


Design 

Requirement 

Design  Goal 

Lamellar  Rhomb 

Capabilities 

(PRR) 

Lamellar  Rhomb 

Capabilities 

(PCR) 

Sinusoidal  Rhomb 

Capabilities 

(PCR) 

Rhomb  efficiency 
for  incoupled 
oscillator  beam 

Greater  than  90% 

(or  sufficient  to 
saturate  the 
amplifier  output) 

81% 

79% 

38%  (must  saturate 
amplifier  output) 

Rhomb  efficiency 
for  outcoupled 
high-power  beam 

Greater  than  95% 

81% 

92% 

96% 

Rhomb  efficiency 
for  orthogonal 
polarization 
return  beam 

Less  than  5% 

0.04% 

0.04% 

0.09% 

effort  in  ion  milling  deep  grooves  (1-micron  depth)  indicates  that 
there  are  materials  control  problems  (surface  crazing,  whiskering, 
uneven  milling)  which  need  to  be  overcome. 

2.  The  grating  efficiency  curves  for  a  sinusoidal  shallow  groove 
design  vary  slowly  with  groove  depth  and  incident  angle.  This 
translates  into  less  stringent  demands  on  fabrication  of  groove 
depth  uniformity  and  system  angle  tuning.  A  sinusoidal  test 
grating  will  be  procured  to  verify  efficiency  performance  over  the 
spectral  band  of  interest. 

5.4  QUARTER-WAVE  RETARDER  STUDY 

The  objective  of  this  effort  was  to  evaluate  coatings  for  the  P](7), 

Pj  (8) ,  P 2(7)  and  P£(8)  laser  lines  that  yield  single-pass  phase  retardance 
of  90  degrees.  The  approach  investigated  was  a  double-bounce  reflection 
per  pass  (Figure  5-34).  The  analyses  studied  were  results  from  the  McLeod 
Code  (University  of  Arizona),  and  a  propriety  code  of  OCLI.  The 
overarching  physical  requirements  throughout  this  investigation  were  the 
use  of  phase  retarding  materials  of  proven  high  optical  fluence  damage, 
high  reflectance  at  the  wavelengths  of  interest,  high  resistance  to 
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Figure  5-34.  90-deg  retardance 
coating  system. 

abrasion,  good  adhesion  characteristics,  and  good  process  control abi 1 i ty 
tor  uniform  deposition,  nominally  1  percent  random  thickness  error  per 
layer. 

It  should  be  noted  that  work  is  ongoing  for  0-deg  phase  retardance 
coatings  for  ALPHA  II  under  the  NWC  HF  coatings  program. 

5.4.1  Pi scussion 

Figure  5-35  shows  a  typical  curve  of  phase  shift  versus  wavelength  for 
an  incident  angle  of  45  degrees  using  a  multilayer  approach.  "L"  repre¬ 
sents  a  relatively  low  refractive  index  materials  such  as  sapphire;  "H" 
represents  a  high  refractive  index  material  such  as  zinc  sulfide;  and  "m" 
represents  the  number  of  layers  of  films,  usually  number  4-5.  Both  s  and  p 
polarizations  have  their  own  characteristics  signatures.  The  high  disper¬ 
sion  of  the  curves  make  this  unacceptable  for  broadband  applications. 

Figure  5-36  shows  the  effect  of  changing  the  outermost  layer  to  a  lew 
refractive  index  material  on  the  phase  retardance  curve:  the  slope  of  the 
dispersion  turns  negative.  When  used  together  with  a  "positive"  design 
such  as  in  a  double-bounce  configuration,  the  net  dispersion  over  the 
spectral  band  of  interest  can  be  significantly  decreased. 

Phase  retardance  curves  for  net  0-degree  phase  as  a  function  of  angle 
of  incidence  are  shown  in  Figure  5-3/.  Curve  A  (angle  of  incidence  =  50 
degrees)  is  equivalent  to  what  might  be  expected  of  a  90-degree  phase 
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retardance  curve.  Thus,  it  is  estimated  that  the  bounce  retardance  error 
for  a  double-bounce  single-pass  retarder  set  would  be  about  6  degrees,  or  a 
round-trip  retardance  error  of  12  degrees;  that  is,  the  return  beam  to  the 
beam  polarization  isolator  would  be  slightly  elliptical  with  a  power  compo¬ 
nent  in  the  orthogonal  direction  of  about  1  percent  back-propagati ng  to  the 
oscillator  subsystem. 

For  the  material  currently  used  for  the  films,  the  mechanical  coating 

I 

resistance  was  rated  to  be  moderate,  and  adhesion  characteristics  would 
conform  to  MIL-C-48497A  (slow  tape  pull).  The  damage  threshold  for  power 
fluence  is  TBD,  however  the  estimate  by  material  type  indicates  that  it 
would  be  on  the  order  of  ALPHA  1  damage  threshold.  Ihe  average  reflectaiue 
was  estimated  to  be  0.992. 
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5.4.2  Summary 


Practical  90-degree  phase-retardance  coating  systems  can  be  made  with 
small  dispersion.  Damage  threshold  needs  to  be  experimentally  verified  for 
high-power  systems.  Design  goals  for  reflectance  can  be  optimized  for 
specific  geometries  and  systems,  i.e.,  angle  of  incidence  and  bandwidth. 
Coating  designs  are  applicable  to  most  HtXs  by  standard  evaporation  methods 
which  will  make  the  component  lower  risk. 


6.  BEAM  DIRECTOR  SUBSYSTEM 


6.1  BEAM  DIRECTOR  PRELIMINARY  CONCEPT  AND  TRADE  ISSUES 

6.1.1  Objective  and  Approach 

The  objective  of  the  beam  director  design  effort  is  to  develop  the 
analytical  tools  ^nd  approaches  with  which  to  generate  a  preliminary 
configuration  that  satisfies  all  the  mission  requirements  and  guidelines 
with  minimum  design  complexity. 

The  approach  was  to  first  use  geometrical  ray  traces  for  the  initial 
beam  director  configuration,  which  would  result  in  several  candidate  design 
approaches.  Then  the  CODE  V  computer  program  was  to  be  used  to  more  accu¬ 
rately  model  and  calculate  the  wavefront  error  induced  by  the  candidate 
design  approaches.  This  would  allow  selection  of  the  design  approach  with 
lowest  overall  OPD  for  the  required  small  angle  pointing. 

Another  activity  is  to  address  the  piston  and  tilt  tolerances  on  the 
primary  mirror  segments  by  using  a  computer  code  to  propagate  the  primary 
mirror  piston  and  angular  errors  through  the  amplifier  and  SBS  cell  to 
determine  allowable  segment  misalignment.  Results  from  both  analyses  are 
presented  in  the  following  sections  of  this  report. 

The  Beam  Director  Subsystem  (BDS)  represents  a  challenging  design  task 
because  of  the  conflicting  requirements  for  this  subsystem.  The  BDS  must 
provide  the  beam  to  targets  within  a  single  field  of  view  with  rapid  retar¬ 
getting  capabilities,  and  yet  must  contribute  small  OPDs  in  the  necessary 
off-axis  pointing  modes.  Also,  the  BDS  contribution  to  overall  system 
jitter  must  be  within  budgeted  requirements,  and  the  beam  obscuration  must 
not  exceed  the  required  value.  Also,  the  requirement  affecting  almost  all 
aspects  of  the  BDS  design  is  the  minimum  brightness  required  of  the  overall 
system,  which  in  addition  to  the  issues  just  mentioned  also  affects  aper¬ 
ture  diameter,  wavefront  error,  and  mirror  power  handling  requirements. 

6.1.2  Beam  Director  Major  Trade  Issues 

Several  areas  have  been  identified  that  require  trade  studies  of 
candidate  design  approaches  to  meet  the  requirements.  These  issues  are 
presented  below. 


6. 1.2.1  Trade  Issue:  Primary  Mirror  f/Number 


Relevant  issues  and  comments:  Separation  distance  between  the 
secondary  and  primary  mirrors  affects  structural  and  optical  considera¬ 
tions.  A  longer  secondary  support  structure  is  less  stable  than  a  shorter 
one,  but  primary  mirror  segment  error  requirements  become  easier  at  greater 
f/numbers.  The  baseline  is  1.0-1.25. 

6 . 1 . 2 . 2  Trade  Issue:  Beacon  Separation  Distance 

Relevant  issues  and  comments:  The  separation  distance  between  the 
free-flying  beacon  mirror  and  the  primary  mirror  affects  both  optical  and 
pointing  considerations.  Larger  separation  distances  result  in  easier 
primary  mirror  segment  alignment  tolerances,  but  the  pointing  accuracy  and 
station  keeping  requirments  are  relaxed  as  the  separation  distance  is 
decreased. 

6. 1.2. 3  Trade  Issue:  Small-Angle  Pointing  Optics 

Relevant  issues  and  comments:  Using  a  spherical  primary  and 
concentric  (zero  field)  optics  allows  small-angle  pointing  with  minimum  OPD 
and  obscuration.  Actual  pointing  is  accomplished  by  moving  the  secondary 
and  tertiary  mirrors  around  a  common  center  of  curvature,  and  removing  the 
resulting  tilt  and  decentration  with  two  movable  relay  optics.  The  alter¬ 
native  is  a  wide  field  of  view  system  that  has  significant  off-axis 
pointing  errors. 

6. 1.2. 4  Trade  Issue:  Single  or  Dual  Secondaries 

Relevant  issues  and  comments:  The  high  and  low  power  beams  each  must 
reflect  off  a  secondary  mirror,  and  this  may  be  a  common  mirror  or  separate 
mirrors.  A  single  secondary  reduces  obscuration  and  allows  conjugation  of 
secondary  mirror  jitter.  Dual  secondaries  produce  more  obscuration  and  do 
not  allow  conjugation  of  jitter,  but  may  yield  more  optical  design  flexi¬ 
bility.  The  baseline  is  a  single  secondary. 

6. 1.2. 5  Trade  Issue:  Muitiseqment  or  Membrane  Primary  Mirror 

Relevant  issues  and  comments:  A  membrane  mirror  composed  of  light¬ 
weight  material  would  offer  weight  and  scaling  advantages,  but  has  not  yet 
been  proven  feasible.  A  multisegmented  glass  mirror  is  the  baseline 
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design,  using  identical  segments  no  larger  than  the  anticipated  state-of- 
the-art  manufacturing  capability. 

6. 1.2.6  Trade  Issue:  Phase  Conjugation  of  the  Primary  Mirror 

Relevant  issues  and  comments:  Tolerances  on  primary  mirror  alignment 
error  are  reduced  by  two  orders  of  magnitude  through  use  of  phase  conjuga¬ 
tion.  Therefore  this  is  the  baseline  approach,  which  has  the  added  advan¬ 
tage  of  a  beacon  mirror  that  drastically  improves  the  beam  quality  of  the 
oscillator  beam  without  spatial  filters,  since  it  must  propagate  a  rela¬ 
tively  long  distance  to  reach  the  beacon.  Figure  6-1  shows  the  significant 
relaxation  of  primary  mirror  segment  alignment  requirements  as  a  result  of 
phase  conjugation  as  a  function  of  beacon  distance  and  f/number. 


RELAXATION  FACTOR  IN  ALLOWABLE 
PHASE  ERROR  AT  PRIMARY  MIRROR 


f  CONJ 
f UNCORR 


PRIMARY 

t/NUMBER 


£Conj  =  allowable  primary  mirror  phase  error 
with  phase  conjugation 

£uncorr  =  allowable  primary  mirror  phase 
error  without  phase  conjugation 

Allowable  error  increases  with  greater  beacon 
distance 

Allowable  error  increases  with  primary 
f/number 

Results  based  on  analysis  of  path  length 
mismatch  due  to  beacon  and  mirror  segment 
piston  error 


L/D  BEACON  DISTANCE 


Figure  6-1.  Conjugated  primary  yields  significant  relaxation  of 
allowable  primary  segment  errors 
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6.1.3  Beam  Director  Preliminary  Configuration  Summar 


The  beam  director  preliminary  configuration.  Figure  6-2,  shows  the 
schematic  side  view,  the  front  view,  and  a  view  of  the  beacon.  Initial 
CODE  V  analysis  has  shown  that  in  order  to  keep  off-axis  pointing  aber¬ 
rations  low,  the  magnification  of  the  high  power  beam  should  be  less  than 
or  equal  to  10.  The  magnification  is  defined  as  the  ratio  of  the  output 
nigh  power  beam  diameter  to  the  diameter  of  the  high  power  beam  at  the 
secondary  mirror.  Also,  the  f/number  of  the  primary  mirror  should  be 
greater  than  or  equal  to  1.0  to  reduce  off-axis  pointing  aberrations; 
beacon-to-primary  separation  distance  of  20  times  the  primary  diameter 
shows  reasonable  OPDs. 

6.2  GEOMETRIC  OPTICS  MULTISEGMENT  BEAM  DIRECTOR  SUBSYSTEM 
6.2.1  Introduction  to  MICROCODE 


The  design,  analysis,  and  tolerancing  of  phased  array  systems  require 
the  ability  to  accurately  model  large  segmented  optical  surfaces  composed 
of  elements  that  may  be  individually  tilted  and  decentered. 

The  primary  lens  design  tool  currently  used  at  TRW  is  CODE  V,  a 
software  package  leased  from  Optical  Research  Associates.  CODE  V  has  the 
capability  to  analyze  only  systems  containing  monolithic  continuous 
surfaces.  It  is  a  leased,  fully  protected  program  which  cannot  be  modified 
by  the  individual  user  for  unique  situations.  Because  of  this  constraint, 
it  was  determined  that  TRW  would  have  to  develop  its  own  analysis 
capability  to  accurately  model  optical  systems  containing  segmented 
surfaces. 

This  software  package  is  MICROCODE,  a  Modifiable  Jnteracti Ve  £ode  for 
Research  in  Optics  by  £omputer  Optical  Design  and  Evaluation.  Its  purpose 
was  to  provide  an  accurate,  modifiable,  analysis  tool  that  could  grow  as 
needed.  Currently  MICROCODE  consists  of  five  main  sections:  (1)  main 
program,  (2)  ray  tracing,  (3)  analysis,  (4)  segmentation,  and  (5)  wavefront 
calculation,  which  are  described  below. 


The  syntax  of  MICROCODE  has  been  kept  as  consistent  as  possible  with 
CODE  V.  MICROCODE  consists  of  a  main  program  with  many  subroutines  that 
carry  out  specific  tasks.  MICROCODE  is  a  bottom-heavy  code,  that  is 
complexity  builds  further  down  the  program.  Essentially  top  subroutines 
call  more  complex  subordinate  subroutines.  The  main  program  includes: 

1.  Data  base  and  its  management.  There  are  two  main  data  arrays: 

LENS  and  IPTR.  LENS  is  the  information  array.  While  reading  in 
surface  information  the  indirect  address  of  the  value  is  saved  in 
IPTR,  with  the  actual  value  being  stored  in  LENS.  When  working  in 
optical  design,  information  has  to  be  readily  accessed  and  modi¬ 
fied.  Indirect  addressing  was  chosen  because  of  the  additional 
flexibility  it  provides. 

2.  Parser.  This  interprets  user  input  for  the  program's  use.  The 
input  is  parsed  for  character  and  numerical  information. 

MICROCODE  can  read  information  from  terminal  or  disk. 


3.  System  Setup.  In  this  area  of  the  program  the  optical  system  is 
set  up.  Refractive  indices  are  determined  for  each  surface  of  the 
system.  MICROCODE  has  the  capability  of  accepting  either  known 
glasses  or  special  glass  runs.  This  allows  the  user  to  exactly 
identify  the  refractive  index  in  question.  Also,  if  exact 
refractive  index  information  is  not  available  for  a  specific  wave¬ 
length,  MICROCODE  can  interpolate  the  given  data  to  get  the 
unknown  index  information. 

In  setup  the  unknown  parameters  of  the  system  are  resolved. 

Any  solves  or  pickups  are  numerically  resolved  and  the  determined 
information  is  placed  in  the  LENS  data  array.  Solves  allowed  in 
MICROCODE  include  marginal  and  chief  ray  angle  solves;  marginal 
and  chief  ray  angle  of  incidence  solves;  marginal  and  chief  ray 
aplanatic  solves,  paraxial  chief  and  axial  ray  height  solves; 
paraxial  marginal  ray  solves  at  known  thickness;  and  center  of 
curvature  solves.  First  order  characteristics  of  the  system  are 
also  determined,  through  a  paraxial  ray  trace. 

6. 2. 1.2  Real  Raytracing 

MICROCODE'S  real  ray  trace  follows  a  general  skew  ray  through  the 
system.  MICROCODE  has  the  capability  of  ray  tracing  through  the  following 
types  of  surfaces:  spherical,  conic,  toric  (cyl i ndi cal ) ,  aspheric, 
anamorphic  aspheres  (nonrotational ly  symmetric),  and  phase  conjugation 
assembl ies. 


A.  BEAM  DIRECTOR  -  SIDE  VIEW 


B.  BEAM  DIRECTOR  -  FRONT  VIEW 
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C.  BEAM  DIRECTOR  BEACON  MIRROR 


Figure  6-2.  Ream  director  concept 


MICROCODE  can  also  handle  individual  surface  tilts  and  decenters. 

Each  surface  can  have  any  combination  of  decenters  in  the  x  direction,  x; 
decenters  in  the  y  direction,  y;  yz  plane  tilts,  a;  xz  plane  tilts,  6;  and 
xy  plane  tilts,  y. 

Normally  MICROCODE  tilts  and  decenters  in  the  following  order 
(x,y ,cx,8,y) ,  following  CODE  V's  convention.  The  user  can  determine  how  the 
tilts  and  decenters  happen  by  specifying  whether  or  not  the  coordinate 
system  of  the  ray  trace  should  change  relative  to  the  surface,  and  in  what 
order  the  tilts  and  decenters  should  occur.  Allowed  forms  of  decentering 
are 

1.  The  coordinate  break  to  start  a  new  coordinate  system 

(DECENTER) 

2.  The  local  coordinate  breaks  for  one  surface  only 

(DECENTER  and  RETURN) 

3.  The  coordinate  break  to  start  a  new  coordinate  system  with 
operations  in  reverse  order  and  sign  (-y,  -8,  -a,  -y,  -x) 

(Ref.  6.1).  (REVERSE  DECENTER) 

6. 2. 1.3  Analysi s 

The  analysis  section  controls  the  geometrical  performance  evaluation 
of  an  optical  system.  The  analysis  section  can  give  any  or  all  of  the 
following  information:  first  order  system  properties,  single  ray  traces, 
ray  traces  of  “fans"  of  rays,  paraxial  aberration  data,  and  surface-by¬ 
surface  printouts  of  real  ray  coordinates  and  angles. 

6. 2. 1.4  Segmentation 

At  this  time  the  user  can  define  up  to  two  segmented  surfaces  per 
optical  system.  Each  of  these  segmented  surfaces  can  have  up  to  100 
individually  defined  segments.  Each  segment  is  defined  by  its  radius,  x 
and  y  coordinates,  and  its  tilts  and  decenters.  Tilts  and  decenters  are 
defined  the  same  as  in  the  main  real  ray  tracing  subroutines,  except  now 
there  is  the  freedom  to  include  piston  errors,  or  z  direction  decenters. 
Editing  of  the  segment  information  can  occur  interactively  as  the  program 
runs . 


6.1  CODE  V  Manual ,  Version  695,  Data-35. 
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The  segmented  surface  ray  trace  controls  the  matrix  transformations  *./ 

necessary  to  accurately  ray  trace  through  a  segmented  surface  and  is  fully 
integrated  into  the  real  ray  tracing  algorithm.  The  procedure  to  ray  trace 
through  a  segmented  surface  includes  identifying  with  which  segment  a  ray 
will  intersect,  applying  the  necessary  coordinate  transfors  to  perform  that 
element's  tilts  and  decenters,  ray  tracing  through  the  element,  and, 
finally,  reestablishing  the  coordinate  system  of  the  rest  of  the  optical 
train. 

6 . 2 . 1 . 5  Wavefront  Calculation 

The  optical  path  difference  algorithm  is  based  on  work  originally  done 
by  H.  H.  Hopkins  in  1952,  and  later  modified  by  Josef  Meiron.  This  OPD 
calculation  is  fully  integrated  into  the  ray  trace  routines,  and  is  done 
automatically  for  every  ray  traced.  At  this  time  MICROCODE  only  performs 
single  ray  OPD  cal cul ati ons . 

6. 2. 1.6  Future  Plans 

The  planned  expansion  for  MICROCODE  includes  adding  further  wave 
optical  evaluation  tools:  OPD  wavefront  mapping,  rms  wavefront  error  ^ 

calculations,  determination  of  the  Strehl  ratio,  calculation  of  the 
diffraction  based  image  of  a  point  object  (point  spread  function),  and 
calculation  of  the  frequency  response  of  an  optical  system  (MTF). 

6. 2. 1.7  Summary 

When  complete,  MICROCODE  will  enable  accurate  prediction  of  the 
performance  of  systems  containing  segmented  surfaces  and  phase  conjugation 
assembl ies. 

6. 2. 1.8  System  Modeled 

The  system  modeled  consisted  of  an  elliptical  primary  and  two  incepen- 
dent  aspheric  secondary  mirrors.  There  was  a  12  to  1  ratio  of  the  diam¬ 
eters  of  the  secondaries  and  the  primary  mirrors.  The  wavelength  used  was 

2.8  microns. 

Light  fills  the  elliptical  primary  from  a  point  source  located  12 
times  the  diameter  of  the  primary  away,  it  is  then  reflected  down  to  an 
aspheric  low-power  secondary  mirror,  and  enters  the  phase  conjugation 
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assembly.  After  exiting  the  phase  conjugation  assembly,  the  light  impinges 
on  the  high  power  secondary  where  it  expands  to  once  again  fill  the  primary 
mirror,  and  then  exits  the  system  collimated. 

6. 2. 1.9  System  Results 

Phase  conjugation  improves  the  system's  performance.  However,  when 
the  same  point  is  not  struck  twice  on  the  primary  mirror,  degradation  in 
performance  occurs,  simply  because  piston  error  has  not  been  fully  can¬ 
celled  by  phase  conjugation.  The  results  support  and  reinforce  results 
found  through  the  use  of  IPAGOS  and  CODE  V.  Phase  conjugation  will  improve 
the  performance  of  the  systems  under  study  but  only  to  the  point  of  ray 
matching  on  the  segmented  primary  mirror. 

Using  MICROCODE  three-dimensional  plots  of  the  OPD,  Figures  6-3  and 
6-4,  were  created  for  the  case  in  which  one  segment  has  a  piston  error  of 
one  wavelength.  Figure  6-3  shows  the  plot  with  the  phase  conjugation 
mirror  in  place,  while  Figure  6-4  shows  the  plot  with  a  plane  reflecting 
mirror  instead. 

6.3  WAVE  OPTICS  ANALYSES  OF  MULTISEGMENT  BEAM  DIRECTOR  SUBSYSTEM- 

MULTISEGMENT  BEAM  DIRECTOR  MISALIGNMENT  ANALYSIS 

This  analysis  was  performed  to  determine  the  sensitivity  of  beam 
quality  of  the  beam  director  output  to  misalignment  of  the  primary  mirror 
segments,  as  a  function  of  beacon  location  and  telescope  f/number,  using  a 
three-dimensional  wave  optics  calculation.  The  sources  of  beam  quality 
degradation  are  assumed  to  be  solely  from  the  beam  director  telescope 
itself  and  optical  path  differences  between  the  low-  and  high-power  beams. 
In  addition,  a  two-dimensional  wave  optics  code  was  to  be  used  to  provide  a 
simulation  of  actual  fields  at  the  entrances  to  the  amplifiers  and  phase 
conjugator.  These  fields,  used  in  conjunction  with  other  codes  which  model 
amplifier  and  SBS  cell  operation,  determine  additional  limits  on  segment 
misalignments  based  on  conjugation  fidelity. 

6.3.1  The  Model 

The  source  of  the  field  into  the  telescope  is  assumed  to  be  c.  point 
source  (beacon)  located  at  a  variable  distance  from  the  primary  and  along 
its  axis.  The  field  from  the  beacon  is  intercepted  by  the  segmented 
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primary,  focused,  and  then  is  collimated  by  the  low-power  secondary  and 
passed  into  the  optical  system  behind  the  primary  mirror.  The  entire 
optical  train  from  the  secondary  to  the  phase  conjugator  is  assumed,  in 
this  model,  to  add  no  aberrations  to  the  field  and  to  contain  no  limiting 
apertures.  Wave  optics  propagation  through  all  the  optical  elements  in  the 
train  can  then  be  described  by  a  single  equivalent  propagation  step  through 
an  effective  distance  dependent  on  the  magnifications  and  separation  of  the 
elements  in  the  train.  The  phase  conjugator  itself  is  assumed  to  be 
perfect  with  no  clipping  of  the  field.  A  second  optical  train  -  not 
necessarily  identical  to  the  first  -  is  traversed  from  the  phase  conjugator 
to  the  high  power  secondary  where  the  beam  is  remagnified  to  fill  the 
misaligned  primary  and  produce  a  collimated  output  beam.  At  this  point  the 
beam  quality  is  evaluated  by  focusing  the  output  beam  and  integrating  the 
intensity  over  a  square  bucket  of  width  equal  to  the  distance  between  the 
first  nulls  at  the  focus  of  a  uniform  field  filling  a  perfectly  aligned 
primary  mirror.  The  ratio  of  the  power  in  the  bucket  of  the  ideal  system 
to  that  of  the  misaligned  system  is  the  beam  quality  number. 

The  primary  and  secondary  mirrors  are  assumed  to  be  circular  para¬ 
boloids.  While  the  actual  mirrors  will  probably  not  be  exactly  para¬ 
boloidal,  the  differences  are  not  expected  to  be  of  great  importance  in 
affecting  segment  misalignment  sensitivity.  Using  the  Fresnel -Ki rchhoff 
integral  representation  of  the  physical  optics  field,  together  with  the 
preceding  assumptions,  allows  the  (high-power)  field  leaving  the  primary  to 
be  evaluated  analytically.  Were  this  not  the  case,  numerical  integrations 
would  be  required  to  propagate  the  field  through  the  optical  system  and, 
considering  the  very  large  Fresnel  numbers  involved  in  the  propagation 
steps,  would  result  in  a  problem  much  more  difficult  to  solve  to  the 
required  degree  of  accuracy.  (A  numerical  propagation  method  would, 
however,  have  the  advantage  of  greater  flexibility  -  for  example,  in  per¬ 
mitting  arbitrary  aberrations  to  be  introduced  within  the  optical  system.) 
Mathematical  details  of  the  model  and  the  derivation  of  the  equations  may 
be  found  in  Section  6. 3. 1.1. 

The  primary  mirror  is  represented  by  an  array  of  contiguous  segments 
whose  projections  in  the  x-y  plane  (i.e.,  the  plane  normal  to  the  optic 
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axis)  are  squares.  The  array  is  truncated  at  the  corners  so  as  to  have  an 
approximately  circular  projection.  The  present  "standard"  configuration  is 
an  array  of  seven  rows  and  columns  with  three  segments  removed  from  each 
corner,  for  a  total  of  37  segments.  The  permitted  misalignments  of  each 
segment  are  any  combination  of  independent  tilts  about  each  of  two  perpen¬ 
dicular  axes  through  its  center  and  a  piston  error,  or  rigid  displacement 
parallel  to  the  normal  to  the  surface.  For  perfectly  aligned  mirror 
segments  and  a  hypothetical  optical  system  for  which  low-  and  high-power 
beams  travel  identical  paths,  the  output  beam  will  be  collimated  and  the 
beam  quality,  by  definition,  will  be  1.  That  is,  the  phase  conjugator  will 
return  a  field  to  each  segment  which  precisely  reverses  all  diffraction 
effects  suffered  on  the  in-bound  path,  and  a  perfect  output  beam  results. 
(In  the  Fraunhofer  zone  each  of  the  segments  has  a  diffraction-limited 
spot  size  of  x/w,  where  w  is  the  width  of  a  segment,  but  with  a  phase 
dependent  on  the  position  of  each  segment.  For  the  perfect  system,  how¬ 
ever,  the  phases  of  the  segment  fields  are  such  that  the  net  field  has  a 
much  smaller  diffraction-limited  spot  size  of  approximately  x/D,  where  D  is 
the  diameter  of  the  primary  mirror.  The  energy  within  the  smaller  spot  is 
used  as  the  measure  of  beam  quality.)  In  general,  however,  the  beacon  will 
be  a  finite  distance  from  the  primary  mirror  and  the  low-  and  high-power 
secondaries  will  be  located  at  different  distances  from  the  primary;  in 
addition,  the  paths  in  and  out  of  the  remaining  optics  may  not  be  quite 
identical.  Thus  the  field  returning  to  the  primary  from  the  phase  con¬ 
jugator  will  not  quite  reverse  the  diffraction  effects  suffered  on  the 
inward  path,  and  the  high-power  field  returning  to  a  particular  segment 
will  experience  some  diffractive  spreading  and  "spillover"  onto  adjacent 
segments.  If  there  are  no  misalignments  of  the  segments,  the  field  lost  to 
one  segment  is  precisely  compensated  by  a  gain  from  the  diffractive 
spillover  onto  it  from  adjacent  segments,  except,  of  course,  for  segments 
at  or  near  the  outer  edge  of  the  primary  which  suffer  a  net  loss.  Accord¬ 
ingly,  the  intensity  at  the  focus  of  the  primary  will  be  reduced  and  the 
beam  quality  number  will  be  somewhat  greater  than  1.0,  even  for  perfectly 
aligned  segments.  The  diffraction  effects  are  characteri zed  by  a  propa¬ 
gation  Fresnel  number  dependent  on  a  net  effective  propagation  distance 
related  to  the  differences  in  optical  paths  of  the  low-  and  high-power 


beams.  Thus  if  the  paths  were  identical,  this  Fresnel  number  is  infinite, 
implying  that  the  diffraction  effects  vanish. 

Suppose  now  that  the  segments  have  random  tilt  misalignments.  The 
phase-conjugated  field  which  returns  to  the  particular  segment  that  it 
encountered  on  its  inbound  path  will  have  a  phase  which  compensates  for  the 
tilt  of  that  segment.  The  diffractive  spillover  onto  adjacent  segments 
will  be  uncorrelated  in  phase  with  the  actual  tilt  of  the  adjacent  segment, 
resulting  in  a  nonzero  net  tilt  of  the  portion  of  the  field  diffracted  onto 
adjacent  elements.  This  field,  when  focused,  then  shifts  away  from  the 
axis  and  the  beam  quality  number  increases.  As  the  segment  tilt  increases, 
the  value  of  BQ  increases  until  all  the  diffracted  spillover  fields  are 
effectively  removed  from  the  central  focal  spot,  at  which  point  the  curve 
of  BQ  versus  tilt  angle  becomes  relatively  flat. 

Another  effect  of  segment  tilt  is  simultaneously  occurring,  however, 
due  to  the  difference  in  optical  paths  of  the  low-  and  high-power  beams. 
Since  the  phase-con jugated  field  of  a  segment  travels  a  different  distance 
than  the  incoming  field,  there  will  be  a  net  displacement  of  the  beam. 

Once  again,  a  spillover  of  the  field  to  neighboring  segments  occurs  and 
that  portion  of  the  field  is  lost  to  the  focal  spot  in  the  far  field.  In 
contrast  to  the  diffraction  effect,  the  beam  quality  degradation  caused  by 
the  geometric  shifting  of  the  beams  continues  to  increase  with  increasing 
segment  tilt  angle. 

Quite  similar  behavior  occurs  for  the  case  of  segment  piston  error. 

The  portion  of  the  field  leaving  a  segment  that  becomes  di ff racti vely 
coupled  to  adjacent  segments  on  its  return  path  has  a  conjugated  phase 
which  is  uncorrelated  with  the  piston  error  of  the  segment  it  actually 
encounters.  The  resulting  phase  aberration  will  quickly  remove  the 
spillover  field  from  the  central  focal  spot  as  segment  piston  error 
increases  from  zero.  In  addition,  even  the  phase  conjugated  field  which 
returns  to  its  segment  of  origin  will  incompletely  compensate  for  the 
segment  piston  error.  This  is  because  the  phase  change  in  the  field 
reflected  from  a  segment  displaced  in  the  direction  of  its  normal  depends 
on  the  direction  of  the  ray  reaching  the  surface.  As  long  as  the  beacon  is 
at  a  finite  distance  from  the  mirror  the  low-power  input  field  and  the 
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phase-conjugated  output  field  will  reach  a  point  on  the  segment  at  dif¬ 
ferent  angles,  with  the  result  that  a  net  uncompensated  piston  error  and  a 
consequent  degradation  in  beam  quality  will  remain.  Consideration  of  the 
relative  positions  of  the  beacon  and  the  low-  and  high-power  secondaries 
suggests  that  beam  quality  will  worsen  for  shorter  beacon  distance  from  the 
primary  and  for  lower  f/numbers.  In  contrast,  the  diffractive  spillover 
effect  on  beam  quality  degradation  will  diminish  with  lower  f/number. 

6. 3. 1.1  Three-Dimensional  Propagation  Equations 


Equations  are  developed  in  this  section  describing  the  propagation  of 
a  field,  originating  as  a  point  source  at  the  beacon  position,  to  the 
segmented  primary,  the  low-power  secondary,  through  an  optical  train  to  the 
phase  conjugator;  the  return  path  from  the  phase  conjugator  is  through  a 
second  optical  train,  not  necessarily  identical  to  the  first,  to  the  high- 
power  secondary,  the  segmented  primary,  and  then  out  of  the  system  as  a 
collimated  beam.  Figure  6-5  is  a  schematic  diagram  indicating  an  equiva¬ 
lent  lens  train  for  the  round  trip  path.  The  direction  of  propagation  is 
taken  to  be  the  z  direction  and  the  x  and  y  axes  lie  in  a  plane  transverse 
to  z.  Each  propagation  step  is  described  by  the  Fresnel  approximation  to 
the  Fresnel -Ki  rchhoff  equation.  The  operations  to  be  performed  at  each 
step  are  enumerated  below. 

Step  1.  Field  at  primary  mirror  due  to  point  source  a  distance  LQ 
from  the  primary 
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Figure  6-5.  Equivalent  lens  train. 


where  EQ  is  an  arbitrary  constant  amplitude  of  the  source,  <5(x,y)  is 
the  Dirac  delta  function,  k  =  2tt/ X ,  and  j  =  /-l.  Thus  the  field  at  the 
primary  is 

JkE0  (  k  2  2  ) 

EpM(x’y) =  exp  j-J  (x  + 

Note  that  the  sign  convention  being  used  corresponds  to  an  exp(+jwt) 
time  dependence  to  the  field. 

Step  2.  Apply  the  primary  mirror  aberration  function  S(x,y).  The 
function  S  represents  the  OPD  variations  due  to  the  misalignments  (tilt, 
piston  error)  of  the  primary  mirror  segments.  This  function  multiplies  Epp| 

l 

EpM(x,>  )  ~  Epjv)(x,y)  S(x ,y ) 

The  explicit  form  of  S  wi 1 1  be  described  later. 

3.  Reflect  off  the  primary  mirror.  The  primary  mirror  is  assumed, 
for  purposes  of  the  misalignment  sensitivity  study,  to  be  a  portion  of 
a  circular  paraboloid.  Thus  the  fiel  d  EPM>  when  reflected  by  the  mirror, 
is  multiplied  by  a  phase  factor  to  give  E pM : 


where  Fp  is  the  focal  length  of  the  primary. 

Step  4.  Propagate  a  distance  to  the  low-power  secondary 

Es(x,y)  =  2^  //  epm(x  1  >y ' )  exP  jiq-  [(x-x')2  +  (y-y')2]{dx,dy' 

where  the  integration  extends  over  the  area  of  the  primary  projected  onto 
the  x,y  plane. 
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Step  5.  Reflect  from  the  low-power  secondary. 

Es(x,y)  =  Es(x,y)  exp  j  ^ —  (xV)  -JkOsLp  x  +  eSLP  y)j 

where  the  secondary  (assumed  to  be  paraboloidal  with  focal  length  F^p)  is 
allowed  to  have  small  tilts  of  e^p,  eSLP  about  the  Y  and  x  axes,  respec¬ 
tively.  FSLp  is  to  be  chosen  so  that  the  field  leaving  the  secondary  is 
collimated;  thus  it  is  determined  by  Fp,  L^,  and  L  .  According  to  the  sign 
convention  Fp  >  0,  F$lp  <  0. 

Step  6.  Propagate  through  an  optical  train  to  the  phase  conjugator. 
The  optical  train  may  be  made  up  of  many  elements  whose  total  effect  on  any 
ray  is  assumed  to  be  representable  by  a  ray  matrix  of  the  form 


In  other  words,  the  collimated  input  beam  is  also  collimated  at  the  phase 
conjugator,  and  the  beam  is  magnified  by  My  on  passing  through  the  train. 

The  same  ray  matrix  is  assumed  to  apply  to  both  x  and  y  coordinate  vari¬ 
ables.  Physical  optics  propagation  through  such  an  optical  train  is 
described  by  the  integral 

oo  i  k  M  ) 

EPC(x,y)  =  'rfjr  //  Es(x'*y')  exP)-J  Tip  [(x’-x/mt)2  +  (y'-y/My)2]  dx'dy' 

provided  there  are  no  limiting  apertures  in  the  optical  train.  It  is 
assumed  that  the  system  would  be  designed  with  minimal  beam  aperturing  in 
order  for  the  output  field  to  most  closely  represent  the  conjugate  of  the 
input  field. 


Step  7.  Phase  conjugate.  For  purposes  of  this  model  the  phase  conju 
gator  is  assumed  to  be  perfect,  so  that  its  effect  is  simply 

1  ★ 

Epc(x,y)  =  Epc(x,y). 

Step  8.  Return  to  the  high-power  secondary  through  an  optical  train 
described  by  the  ray  matrix 


A  B 


C  D 


i/mt  bt 

0  l“L 


My  is  the  same  magni f i cati on  as  for  the  optical  train  in  step  (6);  the 
effective  propagation  distance.  By,  may  be  different  from  By.  The  propa¬ 
gation  integral  for  this  step  is 


Es(*,y)  =  ~  r  //  Epc(x,,y')  exp 


j-S  (K-V)2  +  0',-V)2]|dx'd)' 


Step  9.  Reflect  off  high-power  secondary  as  in  step  (5), 


:s(x,y)  =  Es(x,y)  exp / j  (x2+y2)  -Jk(e^HpX  +  e*Hpy) 


where  F^Hp  <  0  and  is  determined  by  the  magnification  of  the  beam  director 
telescope  and  the  condition  that  the  output  beam  be  collimated. 

Step  10.  Propagate  from  the  high-power  secondary  to  the  segmented 
primary.  If  the  two  mirrors  are  separated  in  z  by  a  distance  L^,  then 


This  last  function,  EpM,  represents  the  field  leaving  the  beam  direc¬ 
tor.  Since  the  objective  is  to  evaluate  beam  quality,  this  field  is  to  be 
focused  and  power-i n-the-bucket  in  the  focal  plane  calculated.  This  pro¬ 
cedure  will  be  described  later. 

Now  consider  the  primary  mirror  segment  aberration  function  S(x,y). 
Figure  6-6  is  a  diagram  of  the  projection  onto  the  x-y  plane  of  the  mirror 
segments.  The  mirror  is  made  up  of  an  array  of  segments  in,  for  example, 
seven  rows  and  columns,  with  three  segments  in  each  corner  removed  so  that 
the  entire  array  is  roughly  circular.  (The  actual  number  of  segments  is  an 
input  variable  to  the  code.)  To  simplify  the  calculations  the  segments 
are  assumed  to  have  identical  square  projections  on  the  x-y  plane,  each 
square  having  a  side  of  length  _w.  The  center  of  each  square  has  coordi¬ 
nates  (xQi,  yQi ) ,  i  =  1,  NS,  where  NS  is  the  number  of  segments.  If  the 
mth  segment  suffers  small  tilts  about  the  y  axis  through  its  center  of  eTx 
and  of  0-j-y  about  the  x  axis  through  its  center,  then  the  OPD  change  of  a 
field  reflected  from  the  segment  corresponds  approximately  to  multiplying 
the  field  by 

s^,(x,y)  •  expj-jk2[e"x  (x-xom)  ♦  e"y(y-yom)]j 

where  (x,y)  is  a  point  on  the  mth  segment. 

For  the  case  of  a  piston  error  of  the  mth  segment  the  OPD  change 
suffered  by  a  field  reflected  from  the  segment  depends  on  the  direction  of 
the  ray  onto  (or  leaving)  the  surface  (Figure  6-7).  The  path  length  change 
At  is 

aL  =  2pm  cos  e 
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where  pm  is  the  displacement  of  segment  m  and  8  is  the  angle  between 
the  surface  normal  and  the  ray  onto,  or  leaving,  the  surface.  If  the  ray 
direction  is  specified  by  the  line  joining  the  point  (x,y)  with  a  point  on 
the  z  axis  a  distance  L  in  front  of  the  mirror,  then  for  a  paraboloidal 
mirror  of  focal  length  Fp,  it  may  be  shown  that 

cose  =  1  -  (x2+y2)/8F2 

where 

T  =  F p / ( 1  -  2Fp/L). 

To  account  for  the  phase  change  in  a  field  reflected  from  a  segment  with 
piston  error,  multiply  by  the  factor 

S^(x,y)  =  exp j-jk2pm[l  -  (x2+y2)/8F2]} 

where  (x,y)  is  a  point  on  the  mth  segment.  The  total  aberration  function 
for  the  surface  is 

S(x,y)  *  I  S  (x,y) 
m 

where  the  sum  is  over  all  segments,  and 

Sm(x,y)  ^(x,y ) ’  if  (x»y)  1s  on  the  mth 
segment 

0  ,  otherwise. 

This  function  S(x,y)  is  the  function  that  appears  in  Steps  2  and  11  above, 
although  its  value  differs  in  the  two  cases.  In  Step  2  the  source  point 
for  a  ray  is  at  L  =  L  (the  beacon  position)  while  in  Step  11  the  rays 


leaving  the  primary  are  collimated  so  that  L  =  «.  In  the  former  case  the 
quantity  F  appearing  in  S(x,y)  is  F  =  F  /(1-2F  /LQ)  while  in  the  latter 
case  it  is  F  =  F  .  Note  also  that  the  function  Sm(x,y)  is  separable  in  x 
and  y,  although  S(x,y)  is  not. 

The  operations  as  prescribed  in  Steps  1  through  11  may  now  be  carried 
out.  At  each  step  an  interchange  in  order  of  integration  is  required;  each 
integrand  contains  complex  exponentials  whose  exponents  are  at  most  qua¬ 
dratic  in  the  integration  variables  so  that  the  integrals  may  be  evaluated 
explicitly.  After  considerable  algebra  the  field  leaving  the  primary 
mirror  from  the  £th  segment  is  obtained  in  the  following  form: 

For  (x,y)  on  segment  i  local  coordinates  are  introduced 

relative  to  an  origin  at  the  center  of  the  segment;  i.e.,  x  =  xQ^  +  5  , 
y  =  yQ£  +  n£,  where  -w/2  _<  5  ,  n£  <_  w/2.  The  field  leaving  the  primary 
from  segment  1  is 
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where  all  superscripts  are  indices  and  not  exponents.  Tne  terms  appearing 
in  this  expression  are 
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plus  a  similar  set  of  quantities  with  x  and  £  replaced  by  y  and  n 
everywhere.  In  addition,  C  is  a  constant  amplitude  factor,  NS  is  the 
number  of  segments  in  the  mirror,  M<-  the  magnification  (>1)  of  the  primary 


\  % 


telescope,  Fpg  the  focal  length  of  the  primary,  and  Fp  =  F pQ / ( 1  -  2Fpg/Lg). 
Finally,  the  function  £(Z)  is  the  complex  Fresnel  integral 

S(l)  =  /  expfj  Z  u2}du. 
o 


As  written  above,  the  quantity  A1  is  assumed  to  be  positive.  If  it  is 
negative,  A1  is  replaced  by  [A1!  and  the  function  £  by  its  complex 
conjugate. 

The  quantity  of  primary  interest  for  assessing  the  effects  of  mis¬ 
alignment  is  not  the  output  field  itself  but  the  beam  quality  number.  In 
order  to  calculate  beam  quality  the  output  field  must  be  focused  and  the 
field  in  the  focal  plane  evaluated.  To  within  a  phase  factor,  the  focal 
plane  field  is  the  Fourier  transform  of  the  field  out  of  the  primary. 

Using  the  Shifting  Theorem  this  transform  may  be  written  as  the  sum  of  the 
transforms  of  the  field  leaving  each  segment.  Each  transform  can  be  eval¬ 
uated  analytically.  The  result  for  the  field  in  the  focal  plane  is 
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where  (xp,yF)  is  a  point  in  the  focal  plane  of  a  focusing  lens  with  focal 
length  f,  and 
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with  r  =  'x1  or  'y ‘ .  The  integral  for  G*1  can  be  evaluated,  with  the 
resul t 
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where 
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and  the  remaining  functions  and  constants  have  been  defined  previously.  As 

before,  if  <  0,  (A^1^2  ♦  (-A^)1/2  and  £  +£  •  Identical  equations  are 
6  1 

found  for  G*  (kyp/f)  by  replacing  ' x '  by  'y *  everywhere  in  the  preceding 
equations. 
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Note  that  diffraction  effects  are  represented  by  the  presence  of  the 
Fresnel  integral  functions  whose  behavior  is  determined  by  the  constant  Ai . 
This  constant  (times  0.5}  plays  the  role  of  a  propagation  Fresnel  number; 
i.e.,  for  small  Ai  diffraction  effects  are  large,  while  if  A^  ♦  »  geometric 
optics  applies.  The  value  of  A^  depends  inversely  on  ^  the  differ¬ 

ence  in  effective  propagation  path  lengths  between  the  ingoing  (to  the 
phase  conjugator)  and  outgoing  fields. 

To  calculate  beam  quality,  a  ' Q=1  bucket'  size  is  found  by  finding  the 
first  null  in  the  x  or  y  directions  at  the  focus  of  a  uniform  field  which 
fills  the  primary  mirror.  The  power  in  a  square  bucket  of  this  size,  for 
the  uniform  field,  is  defined  to  be  the  ideal,  or  reference,  power.  The 
power  in  the  same  bucket  for  an  actual  field  is  then  found  by  integrating 
numerically  the  absolute  square  of  the  field  F(xp,yp).  The  ratio  of  the 
reference  to  actual  powers  is  the  beam  quality  number. 

6. 3. 1.2  Two-Dimensional  Model 

A  two-dimensional  wave-optics  model  was  also  formulated  which  is 
similar  to  the  three-dimensional  model,  but  simpler  in  the  sense  that  the 
field  from  the  beacon  is  traced  only  part  way  through  the  optical  system, 
to  the  position  of  the  amplifier  entrance  and  to  the  S8S  cell.  A  two- 
dimensional  formulation  is  used  because  the  codes  which  model  the  amplifer 
(LCFM  code)  and  the  phase  conjugator  (BOUNCE  code)  are  themselves  two- 
dimensional  models.  The  objective  is  to  assess  conjugation  fidelity  under 
conditions  of  primary  mirror  segment  misalignment.  For  a  system  con¬ 
taining,  say,  seven  amplifiers,  the  fields  into  any  one  amplifier  would 
originate  from  five  or  six  primary  mirror  segments.  Thus,  in  a  linear 
dimension,  the  field  from  only  two  or  three  segments  would  form  the  beam 
into  one  amplifier.  Because  of  the  small  number  of  segments  per  beam, 
segment  misalignment  values  are  simply  assigned  as  input  data  to  the  code, 
rather  than  generating  their  values  by  sampling  a  random  distribution  as  is 
done  in  the  three-dimensional  code.  The  equations  for  the  two-dimensional 
wave  optics  model  are  developed  in  the  following  section. 


6.3. 1.3  Two-Dimensional  Propagation  Equations 

Two-dimensional  fields  which  model,  in  the  physical  optics  approxi¬ 
mation,  the  effects  of  beam  director  segment  misalignment  and  passage 
through  an  optical  train  are  required.  The  procedure  for  development  of 
the  equations  follows  closely  Steps  1  through  6  of  Section  6.3. 1.1,  with 
the  exception,  of  course,  that  the  dependence  on  one  of  the  coordinate 
variables  is  dropped.  In  addition,  since  there  are  multiple  amplifiers  the 
field  off  the  primary  mirror  will  be  divided  at  some  point  so  that  the  beam 
into  a  particular  amplifier  will  originate  from  only  a  few  segments.  The 
model  therefore  traces  the  field  from  a  line  source  to  a  few  (typically  two 
or  three)  misalignment  segments  located  at  a  specified  position  on  the 
primary  mirror,  and  then  to  the  secondary  mirror  and  through  an  optical 
train.  Setting  up  the  equations  and  carrying  out  the  required  integrations 
is  straightforward  though  somewhat  tedious.  The  resulting  expression  for 
the  complex  field  is 
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and  the  remaining  quantities,  except  for  bQ,  are  as  defined  in  Section 
6. 3. 1.1.  The  quantity  bQ  is  the  distance  from  the  axis  of  the  primary 
mirror  to  the  center  of  the  group  of  N  segments  whose  effect  is  being 
considered.  The  x  coordinate  of  the  propagated  field  is  measured  with 
respect  to  the  point  in  the  final  plane  which  is  the  geometric  mapping 
through  the  aligned  telescope  and  optical  train  of  the  center  of  the 
segment  group. 

6.3.2  Calculations  and  Results 

A  computer  code  was  written  to  evaluate  equations  derived  for  the 
far-field  intensity  of  the  collimated  high-power  output  beam,  and  beam 
quality  evaluated.  For  purposes  of  the  preliminary  calculations  the 
optical  paths  from  the  primary  to  the  phase  conjugator,  and  from  the  phase 
conjugator  back  to  the  primary,  were  assumed  to  differ  only  in  the  loca¬ 
tions  of  the  low-  and  high-power  secondaries.  The  segment  misalignment 
tilt  angles  in  both  directions  and  segment  piston  error  were  assumed  to  be 
independent  random  variables  with  a  normal  distribution.  Specific  values 
for  the  misalignments  of  each  segment  were  selected  by  specifying  the 
variance  of  each  variable  (tilts  in  the  x  and  y  directions  were  always 
assigned  equal  variance)  and  then  selecting  values  randomly  from  a  zero- 
mean  normal  population  with  the  prescribed  variance.  For  clarity  in 
showing  trends  and  avoiding  excessive  scatter  of  points,  this  procedure  was 
used  for  one  variance  value;  for  other  values  of  the  variance  the  same  set 
of  misalignment  values  was  then  used  but  scaled  in  magnitude  to  provide  the 
new  variance  value.  An  indication  of  the  spread  in  beam  quality  values  for 
different  random  samples  having  the  same  variance  is  shown  in  the  numerical 
results. 


T0.9 
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RMS  tilt  angles  up  to  10  yrad  were  considered.  It  is  assumed  that 
jitter  frequencies  lower  than  about  200  Hz  are  mechanically  compensated; 
higher  frequencies  are  estimated  to  have  a  most  probable  value  of  around 
0.4  yrad.  RMS  piston  errors  up  to  several  hundred  wavelengths  were  con¬ 
sidered.  Two  values  of  beacon  distance  were  used  -  10  and  100  times  the 
primary  mirror  diameter.  In  most  cases  three  values  of  the  primary  mirror 
f/number  were  used  as  paramters  -  f/8,  f/4,  and  f/2.  While  f/1  would  be  of 
interest,  the  model  would  be  of  doubtful  validity  for  this  case,  princi¬ 
pally  because  of  the  paraxial  assumption.  For  the  case  of  a  beacon  dis¬ 
tance  of  10  mirror  diameters,  only  f/4  and  f/2  cases  were  considered,  since 
at  f/8  the  low  power  secondary  would  be  outside  the  beacon  -  a  situation  of 
little  interest. 

The  results  of  the  computations  are  summarized  in  a  set  of  beam 
quality  versus  misalignment  plots,  and  are  shown  in  Figures  6-8  through 
6-11.  They  may  be  used  to  estimate  the  contribution  to  the  beam  quality 
budget  of  segment  misalignment  for  different  values  of  beacon  distance  and 
f/number  parameters.  It  is  expected  that  the  results  may  be  somewhat 
optimistic,  however,  because  of  the  simplifying  assumptions  made  in  the 
model.  For  example,  imperfect  phase  conjugation  caused  by  nonlinear  gain 
effects  in  the  amplifiers,  beam  aperturing  within  the  optical  train,  or  the 
operation  of  the  SBS  cell  itself  would  increase  the  sensitivity  of  beam 
quality  to  beam  director  misalignments. 

For  the  two-dimensional  model,  the  results  to  date  are  only  in  the 
form  of  plots  of  "typical"  fields  which  might  exist  at  the  amplifier  and 
phase  conjugator  positions  assuming  a  particular  set  of  misalignments  and 
optical  train  elements.  Figures  6 - 1 2a  through  6-12f  show  intensity  and 
phase  plots  for  two-dimensional  fields  entering  the  amplifier  and  SBS  cell, 
for  the  segment  misalignments  shown  in  Table  6-1.  Since  the  field  profiles 
depend  on  the  magnifications  and  path  lengths  of  the  optical  train,  the 
results  shown  should  be  taken  as  merely  suggestive;  these  calculations  may 
be  repeated  as  the  actual  optical  elements,  and  their  positions,  become 
better  defined.  It  appears  likely,  however,  that  because  of  the  large 
demagni f i cati ons  required,  the  propagation  Fresnel  numbers  will  tend  to  be 
rather  low.  For  low  Fresnel  numbers  the  spreading  of  the  beam  beyond  its 
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Figure  6-12.  Field  into  amplifier 
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Figure  6-12. 

Field  into  ampl i f ier 

(continued) 

TABLE  6-1.  Optical  Train  Assumed  for  Field  Calculations. 


From 

To 

Distance* 

Magnification 

Beacon 

Beam  director  primary 

10 

— 

Primary 

Beam  director  secondary 

2.3 

1/12.5 

Secondary 

Telescope 

3.3 

1 

Telescope 

— 

1/17.6 

Telescope 

Amplifier  inner  cone 

06 

1 

Inner  Cone 

Outer  cone 

— 

1/3 

Amplifier 

0.67 

1 

Outer  cone 

Inner  cone 

— 

3 

Inner  cone 

SBS 

0.33 

9 

*  Distance  in  units  of  primary  mirror  diameter 


geometric  size  may  be  substantial,  and  difficulties  may  be  encountered  in 
avoiding  beam  aperturing.  The  entrance  aperture  to  the  amplifie-  annular 
gain  region  was  taken  to  be  1.3  cm  (x  =  -0.65  to  +  0.65  in  Figure  6-12f ) , 
and  that  of  the  SBS  cell  to  be  35  cm  (x  =  -17.5  to  x  =  +17.5  in  Figures 
6-13a  through  6 -1 3d ) .  Beam  displacement  due  to  segment  tilt,  and  inter¬ 
ference  effects  from  beam  overlap,  are  also  quite  pronounced.  For  the 
assumed  geometry  at  least,  the  wave  optics  fields  differ  markedly  from  the 
geometric  optics  fields  (e.g.,  uniform  amplitude  and  phase)  normally  used 
as  input  to  the  LFCM  and  BOUNCE  codes. 
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Figure  6-13.  Field  into  SBS 
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Figure  6-13.  Field  into  SBS  (continued). 
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7.  PHASE -CONJUGATION  SUBSYSTEM 


7.1  PHASE  CONJUGATION  SUBSYSTEM  PRELIMINARY  CONCEPT  AND  TRADE  STUDIES 

7.1.1  Objective  and  Approach 

The  objective  of  the  phase  conjugation  subsystem  effort  is  to  develop 
the  analytical  tools  with  which  to  design  the  SBS  cell,  spectral  separation 
grating,  quarter  wave  retarding  coating,  and  focusing  optics. 

The  approach  was  to  use  the  BOUNCE  and  BRIWON  codes  to  characterize 
the  interaction  of  the  following  parameters: 

1.  SBS  cell  input  power 

2.  SBS  cel  1  reflectivity 

3.  Threshold  power  (required  for  10  percent  reflectivity) 

4.  Beam  focusing  geometry  (near-  and  far-field  overlap) 

5.  Beam  combination  (several  beams  in  a  common  focal  volume) 

6.  Beam  propagation  issues  (thermal  blooming,  self-focusing,  etc.) 

7.  Stokes-seed  back-injection  (to  lower  threshold,  but  accomplish 
multibeam  piston-error  correction) 

8.  SBS  conjugation  fidelity  (measure  of  output  beam  quality) 

9.  SBS  medium  properties  (impurities,  temperature,  pressure,  etc.) 

Many  of  these  parameters  were  examined  using  the  previously  mentioned 
codes,  and  the  results  were  used  to  produce  the  preliminary  configuration 
shown  in  this  report.  These  results  are  presented  in  the  remaining  por¬ 
tions  of  this  section. 

7.1.2  Phase  Conjugation  Major  Trade  Issues 

Several  areas  have  been  identified  that  require  trade  studies  of  can¬ 
didate  design  approaches  to  meet  the  requirements.  These  issues  are  pre¬ 
sented  below. 

7. 1.2.1  Trade  Issue:  Type  of  Cell  Medium,  and  Conditions 

Relevant  issues  and  comments:  The  type  of  gas  used  in  the  SBS  cell 
affects  conjugation  fidelity,  reflectivity,  and  threshold  power.  Pressure 


and  temperature  also  affect  these  parameters.  Analytical  results  indicate 
xenon  at  40  atm  and  300  K  has  the  desired  SBS  characteristics.  COS  experi¬ 
ments  will  address  these  issues. 

7. 1.2.2  Trade  Issue:  Cell  Medium  Purity 

Relevant  issues  and  'omments:  The  level  and  type  of  impurities  in  the 
medium  affect  propagation  phenomena  such  as  thermal  blooming  and  self- 
focusing.  COS  experiments  will  anchor  analytical  calculations. 

7. 1.2.3  Trade  Issue:  Flow  Velocity 

Relevant  issues  and  comments:  Flow  velocity  is  relevant  to  thermal 
absorption  issues  and  phonon  rise  time  concerns,  and  is  intimately  asso¬ 
ciated  with  the  flow  direction.  For  transverse  flow,  the  gas  must  move 
fast  enough  to  limit  thermal  absorption,  but  not  fast  enough  to  smear  the 
phonon  grating  which  exists  on  time  scales  on  the  order  of  the  phonon  rise 
time,  since  this  would  reduce  conjugation  fidelity.  Preliminary  calcula¬ 
tions  show  that  a  range  of  20  to  80  meters/sec  is  acceptable. 

7. 1.2.4  Trade  Issue:  Flow  Direction 

Relevant  issues  and  comments:  Flow  in  the  SBS  cell  that  is  either 
transverse  or  antiparallel  to  the  beam  will  either  not  correct  or  correct, 
respectively,  for  the  SBS  frequency  shift  in  the  Stokes  beam.  Analysis 
(see  Section  5.2.3)  shows  no  shift  compensation  is  required  since  there  is 
small  power  loss  in  the  amplifiers  due  to  the  SBS  frequency  shift. 

Therefore  transverse  flow  is  acceptable. 

7. 1.2.5  Trade  Issue:  Number  of  Cells 

Relevant  issues  and  comments:  Each  line  requires  separate  focal 
volumes  in  the  SBS  medium;  this  may  be  accomplished  through  separation  in  a 
common  cell,  or  by  using  individual  cells  for  each  line.  Separate  cells 
appear  easier  to  fabricate,  use  less  gas,  and  eliminate  potential  focal 
volume  interference  effects  that  might  exist  in  a  common  cell. 

7. 1.2.6  Trade  Issue:  SBS  Telescope  f/Number 

Relevant  issues  and  comments:  The  focal  length  divided  by  the  diam¬ 
eter  of  the  focusing  optics  affects  conjugation  fidelity  and  cell  length 
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and  width  requirements.  Probable  range  is  f/10  to  f/50.  COS  experiments 
will  address  this  issue. 

7.1.2  Preliminary  Configuration 

The  four  side-by-side  SBS  cells  were  shown  in  Figure  4-7.  Note  that 
the  fraction  of  energy  that  is  transmitted  through  the  cells  is  directed 
backward  toward  the  rear  of  the  spacecraft  and  out  into  space  so  that  it 
will  not  disturb  any  satellite  hardware.  The  exhausted  helium  and  xenon 
gases  also  are  directed  out  towards  the  back  of  the  spacecraft. 

7.2  PHASE  CONJUGATION  ANALYSIS 
7.2.1  SBS  Modeling 

Phase  conjugation  is  the  wavefront  process  by  which  a  beam  retraces 
its  path  through  aberrations  and  corrects  its  wavefront  distortion  by 
virtue  of  phase  sign  reversal  over  that  of  the  forward  beam.  An  important 
method  of  phase  conjugation  is  SBS  reflection  from  a  phase  conjugation  cell 
containing  an  appropriate  nonlinear  optical  medium.  For  PALS  this  is  xenon 
at  about  40  atm  pressure. 

SBS  phase  conjugation  occurs  by  means  of  the  preferential  amplifica¬ 
tion  of  the  backscattered  wave  that  is  the  phase  conjugate,  i.e.,  the 
(backward-going  Stokes)  wave  with  the  apparent  time-reversed  propagation 
properties  to  those  of  the  incident  wave  (pump  beam). 

In  the  case  of  a  pump  beam  focused  into  the  cell,  the  intense  focal 
region  comprises  a  substantial  part  of  the  so-called  interaction  region  for 
the  SBS  phase  conjugation.  However,  as  seen  in  Figure  7-la,  the  Stokes 
intensity  is  greatest  forward  of  the  focal  region,  and,  in  this  case  at 
twice  threshold  power,  is  restrained  by  the  forward  position  of  the  mathe¬ 
matical  cell  window  at  1  cm  ahead  of  the  focal  center  at  the  origin.  The 
f/number  (lens  focal  length  divided  by  hypergaussi an  beam  diameter)  of  the 
pump  beam  is  40.  The  Rayleigh  range,  i.e.,  the  full  length  of  the 
diffraction  limited  focal  region,  is  approximately 

L  :  4X  ( F / « ) 2  .  1.5  cm  (7.1) 

a  value  which  agrees  with  the  contour  plot  of  Figure  7-lb. 
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Details  of  the  SBS  phase  conjugation  process  are  calculated  in  the 
BOUNCE  code.  The  electromagnetic  wave  equations  for  the  pump  and  Stokes 
beams  are  solved  by  a  combination  of  amplification  and  propagation 
algorithms  working  over  a  discrete  two-dimensional  space.  The  gain  sheets 
at  which  the  amplification  or  gain  is  implemented  are  specified  along  the 
axial  propagation  direction.  Diffractive  propagation  sideways  is  imple¬ 
mented  by  a  fast  Fourier  transform  technique  in  one  transverse  direction. 
Effects  in  the  other  transverse  dimension  are  treated  by  simple  geometrical 
focusing  factors  down  to  a  diffraction-limited  cutoff  distance  at  the  focal 
region. 

The  amplification  derives  from  the  SBS  or  stimulated  Brillouin  scatter 
in  which  pump  photons  are  converted  to  Stokes  photons  at  a  rate  that  is 
proportional  to  the  number  of  photons  of  both  pump  and  Stokes  waves.  The 
spatial  and  temporal  beats  between  the  two  waves  drive  an  acoustic  wave, 
i.e.,  phonons.  This  acoustic  grating  is  the  scatter  mechanism  whereby  the 
complicated  far-field  pattern  of  pump  intensity  spots  is  replicated  by 
Stokes  backscatter.  Aberrations  have  distorted  the  incident  pump  wavefront 
to  produce  this  complicated  far-field  intensity  pattern. 

When  the  phonons  have  a  large  damping  factor,  as  they  usually  do  for 
backward  light  scatter  at  infrared  or  shorter  wavelengths,  the  acoustic 
wave  amplitudes  are  driven  locally  by  the  product  of  the  pump  and  Stokes- 
complex-conjugate  electric  fields.  The  acoustic  wave  frequency  and  wave 
vector,  in  the  plane  wave  decomposition  of  all  the  waves,  are  the  appro¬ 
priate  differences  of  the  pump  and  Stokes  quantities,  and  they  are  also 
related  to  the  acoustic  wave  dispersion  relation  involving  the  speed  of 
sound  in  xenon. 

The  local  nature  of  the  relationship  between  acoustic  wave  amplitude 
and  pump-Stokes  product  is  lost  when  a  cross-flow  exists.  A  convolution  of 
this  electric  field  product  and  the  exponential  function  for  convection  and 
damping  with  transverse  distance  along  the  flow  replaces  this  local  rela¬ 
tionship.  With  increased  cross-flow  speed  the  smearing  effect  of  this 
convolution  is  evident  in  a  decrease  in  gain  coefficient  and  in  a  new  non¬ 
linear  refractive  index  contribution.  The  gain  coefficient  reduction  is  a 
serious  matter  when  competition  with  amplification  for  Stokes  waves  other 


than  the  conjugate  Stokes  wave  is  at  stake.  Typically  the  phase  conjugate 
Stokes  wave  has  about  a  factor  of  two  larger  gain  coefficient  than  for  the 
similar  but  nonconjugate  Stokes  waves,  in  the  absence  of  cross-flow. 

If  the  distance  d  characterizes  the  transverse  distance  over  which  the 
pump-Stokes  product  is  correlated,  then  the  gain  coefficient  is  reduced  by 
the  approximate  factor 


(1  +  2Vr/d) 


(7.2) 


where  r  is  the  reciprocal  damping  rate  constant  for  the  acoustic  wave. 

This  quantity  has  already  been  referred  to  as  the  small-signal  phonon  rise 
time  or  phonon  lifetime.  The  cross-flow  or  transverse  flow  speed  is 
denoted  by  V. 

The  rise  time  for  phonons  generated  in  forward  SBS  is  much  larger  than 
for  the  backscatter  SBS  phonons.  Thus,  the  cross-flow  speed  is  chosen  to 
flush  pnonons  generated  in  forward  SBS  but  not  in  backward  SBS. 

The  BOUNCE  code  is  difficult  to  run  successfully  at  high  reflectivity. 

A  successor  to  BOUNCE  has  been  developed  that  reaches  convergence  in  its 
iterations  in  a  somewhat  different  manner,  and  usually  the  costs  to  run 
BRIWON  (or  Bl)  are  less. 

Instead  of  using  conservation  of  energy  between  the  pump  and  Stokes 
beam  as  a  criterion  of  convergence,  BRIWON  incorporates  energy  conservation 
within  the  amplification  algorithm  at  the  gain  sheets.  For  convergence 
BRIWON  requires  the  same  results  whether  integrating  the  simultaneous  equa¬ 
tions  in  the  forward  direction  or  in  the  backward  direction.  Meanwhile, 
the  pump  boundary  condition  at  the  front  of  the  cell  and  the  Stokes  bound¬ 
ary  condition  of  growth  from  Stokes  noise  mainly  from  the  back  of  the 
interaction  region  are  satisfied. 

Figure  7 - 1 c  also  shows  reflectivity  curves  for  xenon  at  30  and  40  atm. 
The  horizontal  scale  is  input  pump  power  in  units  of  threshold  power  for  the 
40  atm  case.  Threshold  is  defined  as  that  power  resulting  in  10  percent  SBS 
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reflectivity.  The  threshold  power  for  40  atm  is  calculated  to  be  about  150 
kW.  This  is  close  to  that  calculated  from  the  formula  that  follows: 


gIL  =  gPl/(2XF,'#)2  »  25 


(7.3) 


or 

P  =  25X/g 


(7.4) 


where  g  is  the  gain  coefficient,  P  is  the  input  power,  F/fr  is  the  beam 
f/number,  X  is  the  wavelength  and  L  is  the  Rayleigh  range  discussed  above. 

A  round  beam  and  a  somewhat  higher  value  for  the  gIL  product,  namely  30 
instead  of  25,  yield  240  kW,  as  noted  by  NRL's  Whitney,  et  al.  (Ref.  7.1). 

The  horizontal  scale  in  Figure  7 - 1 c  is  logarithmic.  The  two  curves 
for  the  different  pressures  lie  almost  a  fixed  horizontal  distance  apart, 
indicating  that  the  reflectivities  are  in  a  fixed  ratio  to  each  other. 

This  ratio  is  close  to  the  square  of  the  ratio  of  pressures.  This  is 
expected  from  the  scaling  of  the  gain  coefficient,  which  is  proportional  to 
the  square  of  the  pressure.  One  power  of  pressure  is  the  usual  one  from 
the  dependence  on  the  number  density  of  scatters,  while  the  second  power  of 
pressure  is  from  the  phonon  rise  time.  That  is,  the  gain  coefficient  is 
directly  proportional  to  the  phonon  lifetime,  which  is  inversely 
proportional  to  the  phonon  bandwidth.  The  bandwidth  decreases  with 
pressure  because  the  kinematic  viscosity  and  thermal  diffusivity  that  damp 
the  acoustic  waves  for  the  xenon  both  decrease  with  pressure  at  a  given 
temperature. 

The  simultaneous  SB5  phase  conjugation  of  two  beams  has  been  calcu¬ 
lated  with  both  codes.  The  results  are  not  fully  converged.  BR1W0N 
indicates  that  there  is  no  penalty  in  reflectivity  at  a  given  power  with 
two  beams  compared  to  one,  provided  that  there  is  overlap  of  interaction 

volumes  in  the  cel  1 . 

Beam  quality  as  well  as  conjugation  fidelity  characteri sti cs  are  output 
of  the  codes  or  of  ancillary  codes. 


.;•>  7.1  "SBS  with  an  HF  Laser,"  S.  Whitney,  et  al,  Post-deadline  paper, 
presented  at  International  Conference  on  Lasers  1985,  Dec  2-6, 
Las  Vegas,  NV. 
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Nonlinear  processes  that  might  compete  with  SBS  include  stimulated 
Raman  scattering,  nonlinear  absorption,  stimulated  thermal  Brillouin 
scattering,  and  stimulated  Rayleigh  line  scattering. 

A  process  that  competes  with  SBS  is  SRS,  stimulated  Raman  scatter, 
which  is  possible  for  molecular  xenon  formed  at  these  high  pressures  and 
densities.  At  40  atm  about  15  to  20  percent  by  number  is  diatomic  xenon. 
The  conversion  of  pump  to  Stokes  photons  in  SRS  excites  internal  rotation- 
vibration  states  (analogous  to  so-called  optic  phonons). 

SRS  has  a  gain  coefficient  that  is  usually  much  smaller  than  for  SBS. 

A  calculation  based  on  pol ari zabi 1 i ty  data  (Ref.  7.2)  confirms  this.  If 
the  xenon  were  entirely  in  diatomic  form,  i.e.,  as  van  der  Waals  molecules, 
the  SRS  gain  coefficient  at  40  atm  would  be  about  3  x  10-5  cm/MW  compared 
with  the  SBS  gain  coefficient  of  about  5  x  10‘2  cm/MW. 

Competition  with  SBS  might  occur  in  the  short  pulsed  case,  for  which 
the  slower  rise  time  of  the  acoustic  wave  puts  SBS  in  the  transient  regime. 
Then,  the  transient  SBS  gain  coefficient  might  be  less  than  the  gain 
coefficient  for  SRS,  which  might  be  in  the  steady-state  regime  owing  to  the 
much  shorter  rise  time  (about  10~10  s)  for  the  rotation-vibration  modes  of 
the  molecular  xenon. 

Nonlinear  absorption  by  multiphoton  excitation  of  xenon  electronic 
states  is  negligible.  More  than  19  HF  photons  are  needed  to  reach  the 
first  excited  state. 

Impurities  in  xenon  such  as  water  vapor,  carbon  dioxide,  and  nitrous 
oxide  are  to  be  held  to  low  concentrations,  i.e.,  less  than  0.1  ppmv.  The 
linear  absorption  coefficient  is  then  about  3  x  }0'6  cnrl  (that  is,  0.3  per 
km).  The  STBS  or  stimulated  thermal  Brillouin  scatter  process  has  a  gain 
coefficient  (in  the  usual  units  of  cm/MW)  about  700  times  the  absorption 
coefficient  (in  cm'*).  This  yields  a  coefficient  about  5  percent  of  that 
of  SBS,  for  xenon  at  40  atm  and  impurities  at  less  than  0.1  ppmv. 


7.2  M.H.  Proffitt,  et  a  1 .  Can.  J.  Phys.  59,  1459-1480  (1981). 


STBS  results  from  the  absorption  of  energy  from  the  pump  beam,  as  in 
thermal  blooming;  this  energy  causes  a  thermal  volume  expansion  that  drives 
an  acoustic  wave  analogous  to  that  in  SBS.  Characteri sti cs  of  the  Stokes 
backscatter  in  frequency  shift  (107  MHz)  and  phonon  rise  time  (260  ns)  are 
also  similar.  Distinctions  are  seen  only  at  high  frequency  resolution, 
where  STBS  gain  coefficient  maximum  is  offset  from  SBS  maximum  by  about  the 
phonon  damping  ha1 f-bandwidth. 

SRLS  is  stimulated  Rayleigh  line  scatter,  which  is  also  dependent  on 
the  linear  absorption  coefficient.  However,  in  this  process  the  driven 
entropy  "wave"  rather  than  the  driven  acoustic  wave  is  the  scatter  mecha¬ 
nism.  The  gain  coefficient  for  SRLS  is  comparable  to  that  of  STBS,  the 
rise  time  is  somewhat  longer,  however.  Thus,  in  pulsed  cases  SRLS  might  be 
less  important.  However,  in  CW  or  quasi-CW  cases,  SRLS  might  dominate. 

The  frequency  shift  between  pump  and  Stokes  signals  would  be  very  small  for 
SRLS,  compared  to  that  for  STBS,  since  an  acoustic  wave  is  not  involved 
with  SRLS. 

Phase  conjugation  by  backscattered  SBS  also  competes  with  near-forward 
SBS.  The  gain  coefficient  increases  as  the  angle  of  scatter  decreases. 
However,  the  rise  time  for  the  phonons  generated  in  the  near-forward  SBS 
increases  even  more  strongly  with  decreasing  angle  of  scatter.  Thus,  only 
very  long  pulsed  or,  of  course,  CW  cases  are  concerned  with  near-forward  SBS. 

Near  the  conclusion  of  the  PALS  program,  a  time-dependent,  quasi  one¬ 
dimensional  model  has  been  developed  to  set  engineering  limits  on  the  speed 
of  crossflow  or  beam  slewing  for  suppressing  near-forward  SBS  while  not 
diminishing  backward-going  SBS.  The  effective  gain  coefficient  for  various 
angles  of  scatter  can  be  shown,  as  a  function  of  elapsed  time,  to  approach 
the  usual  steady-state  coefficient. 

The  geometry  of  the  focal  region  needs  to  be  incorporated  into  the 
model  so  that  the  SBS  output  at  various  angles  can  be  estimated  as  a  func¬ 
tion  of  time.  The  time  may  be  identified  with  the  ratio  of  the  focal  spot 
radius  to  the  relative  transverse  speed  due  to  crossflow  or  scanning  by  the 
beam. 
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The  well-known  behavior  of  the  scattered  field  according  to  a  Bessel 
function,  whose  argument  is  proportional  to  the  square  r,  *-  of  the  the  prod¬ 
uct  of  the  elapsed  time  and  the  interaction  length,  is  the  basis  of  these 
effective  gain  coefficient  curves. 

7.2.3  Thermal  Blooming,  Self-Focusing,  and  IR  Gas  Breakdown  in  SBS  Cell 

The  effects  on  SBS  phase  conjugation  for  single-pulse  and  quasi-CW 
cases  a>~e  analyzed  by  approximate  wave  propagation  formulae  and  the  BOUNCE 
code  *.ith  modifications.  For  thermal  blooming,  the  quasi-CW  and  pulsed 
case,  are  considered  and  impurity  level  limits  are  set.  In  the  CW  case,  a 
range  of  transverse  flow  speeds  is  specified.  In  the  pulsed  case,  a  limit 
on  time  variation  of  the  beam  deformation  parameter  during  the  phonon  rise 
time  is  set.  For  self-focusing  in  the  quasi-CW  case,  maximum  power  is 
determined  by  the  breakdown  threshold  in  clean  xenon.  In  the  pulsed  case, 
focal  motion  is  limited  to  1/e  of  interaction  length  in  a  phonon  rise  time. 
Combined  thermal  blooming  and  self-focusing  causes  partial  cancellation 
with  an  exponential  intensity  profile. 
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Thermal  blooming  is  the  distortion  of  a  beam  by  the  nonuniform 
decrease  of  refractive  index  produced  by  thermal  volume  expansion  due  to 
absorbed  beam  energy.  Self-focusing  is  the  contraction  and  related 
distortion  of  a  beam  by  nonuniform  refractive  index  increase  due  to  the 
nonlinear  effects  of  the  beam  intensity.  Optical  breakdown  at  HF 
wavelengths  is  the  creation  of  plasma,  usually  in  the  vicinity  of  a  beam 
focus  by  the  acceleration  of  seed  electrons  in  the  beam  electric  field,  the 
impact  ionization  of  the  medium,  electron  avalanche,  and  HF  infrared 
absorption . 


Phase  conjugation  by  SBS  in  the  baseline  xenon  cell  at  40  atm  could 
lose  conjugation  fidelity  and  SBS  reflectivity  if  the  thermal  blooming  of 
the  pump  and  Stokes  beams  were  large  enough.  Rays  of  the  pump  beam  that 
are  bent  away  from  the  SBS  interaction  region  to  where  intensities  are  not 
high  are  not  phase  conjugated.  Thus,  aberrations  introduced  by  thermal 
blooming  cannot  be  wholly  corrected  by  SBS  phase  conjugation. 

The  thermal  blooming  discussed  here  is  caused  by  the  absorption  of 
energy  at  HF  wavelengths  owing  to  impurities  such  as  water  vapor,  carbon 
dioxide,  and  nitrous  oxide.  In  the  CW  or  quasi-CW  cases,  a  transverse  flow 
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is  introduced  to  carry  away  the  SBS  medium  that  is  disturbed  by  the  heating 
effect  of  infrared  absorption.  Without  flow  or  an  equivalent  slewing  of 
the  beam  there  would  be  a  continual  heating  of  the  medium  and  detrimental 
distortion  of  the  beam  until  very  long  times  when  thermal  conduction  would 
set  up  a  steady  state. 

Self-focusing  in  xenon  is  possible  due  to  the  nonlinear  effect  of 
electrostriction  on  the  index  of  refraction.  The  coefficient  that,  charac¬ 
terizes  this  effect  is  large  owing  to  the  large  polarizability  of  xenon 
atoms.  Self-focusing  is  often  accompanied  by  optical  breakdown  in  the 
intense  focal  region.  Breakdown  competes  with,  but  does  not  necessarily 
el imi nate  SBS. 

The  pulsed  cases  of  thermal  blooming  and  self-focusing  are  different 
from  the  CW  or  quasi-CW  cases  because  of  the  shortness  of  the  elapsed  time 
compared  to  the  transit  time  of  sound  waves  over  a  local  beam  radius.  For 
the  focused  geometry  situation,  the  beam  radius  is  taken  to  be  the 
diffraction  limited  spot  radius.  This  size  is  characteristic  of  the  finest 
structure  of  the  strongest  intensities,  for  which  the  thermal  blooming 
effects  are  by  far  the  greatest. 

For  short  pulses,  because  the  medium  cannot  easily  adjust  its  density, 
and  hence,  its  refractive  index,  when  the  sound  waves  have  not  crossed  the 
beam  to  signal  an  adjustment,  the  effects  of  thermal  blooming  and  self- 
focusing  are  smaller  than  for  CW  or  quasi-CW  by  the  square  of  the  ratio  of 
the  elapsed  pulse  time  to  the  sonic  transit  time.  Quasi-CW  refers  to  the 
cases  of  long  enough  laser  pulses  that  the  regions  of  greatest  refractive 
index  change  are  in  their  steady  state  with  regard  to  the  sonic  transit 
time. 

The  effects  on  the  beam  quality  in  either  phenomenon  are  conveniently 
summarized  in  the  distortion  parameter  OP.  This  parameter  accounts  for  the 
change  in  intensity  on  axis  due  to  the  sideward  spread  or  contraction  of 
the  beam.  The  lowest  order  wavefront  distortions  of  piston  and  tilt  are 
correctable  by  phase  conjugation  in  the  usual  manner,  since  changes  in 
relative  intensity  are  not  involved,  though  details  of  location  of  the 
interaction  region  are. 
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The  next  order  wavefront  distortion  of  spherical  focus  or  defocus 
displaces  the  point  of  geometric  focus.  Only  the  time  dependence  of  this 
effect  is  important,  and  only  when  there  is  substantial  change  in  intensity 
distribution  in  the  interaction  region  for  SBS  during  a  so-called  phonon 
lifetime,  which  lifetime  is  the  small-signal  acoustic  wave  rise  or  decay 
time. 

In  the  CW  or  quasi-CW  cases,  the  cross-flow  causes  the  pump  and  Stokes 
beams  to  deflect  against  the  flow  for  thermal  blooming.  The  distortions  in 
the  direction  perpendi cular  to  both  the  flow  and  the  original  beam  direc¬ 
tions  are  essentially  free  of  this  type  of  wavefront  tilt.  For  conve¬ 
nience,  the  distortion  parameter  is  calculated  along  this  perpendicular 
di rect i on . 

It  is  noted  that  thermal  blooming  does  not  merely  defocus  the  beam. 

The  central  region  of  the  beam  intensity  distribution  acts  somewhat  in  this 
manner  while  the  wing  region  acts  oppositely,  in  accordance  with  the 
opposite  sense  of  curvature  of  the  distribution  and  of  that  induced  in  the 
wavefront.  This  is  cvl indrical ly  symmetric  with  no  cross-flow.  In  effect, 
cross-flow  macs  the  beam  intensity  distribution  to  new  positions  stretched 
along  the  flow  direction. 

The  energy  conservation  equation  can  be  used  to  express  the  intensity 
I  on  axis  In  terms  of  the  refractive  index  be  as  follows: 


H  *  Hr  5  I  nds 


I  a  e 


at  ri  =  0 


(7.5) 


where  DP,  the  distortion  parameter,  is  given  by 


n(s,rJ_ ) 


f  ds '  n(s  rj  +  ^n(s  r)  .  ^ln  Ijsr)^ 


(7.6) 


with  the  variables  of  integration  s'  and  s  being  distances  along  ray  paths, 
and  rA  the  transverse  coordinate,  as  the  rays  are  brought  to  a  focus  in  the 
SBS  medium.  The  propagation  distance  in  the  medium  to  the  front  of  the 
focal  region  is  f . 
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The  distortion  parameter  is  easily  evaluated  for  a  Gaussian-shaped 
beam.  The  result  for  the  CW  or  quasi-CW  cases  is  as  follows: 


dA  I6t 
dt  .2 


(7.7) 


where  77  = 
dt 


dn  gP 
dT  n pC 


and  t 


flong-pulse  duration  or 
IXF'/V,  with  cross-flow  V, 


In  the  CW  or  quasi-CW  cases,  the  useful  criterion  is  that  DP  be  less 
than  unity.  The  conclusion  for  thermal  blooming  is  that  the  absorbed  power 
per  unit  cross-flow  speed  be  less  than 


a  PF*/V  <  3  (W/cm)  '  (ni/ s)  (7.8) 

F -  =  beam  F-number 
n  =  refractive  inder 
X  -  wavelength 
V  =  cross-flow  speed 
Cp  =  speci f i c  neat 
P  =  input  power 
a  -  absorption  coefficient 
p  -  density 
T  =  temperature 

The  condition  that  the  distortion  parameter  be  less  than  unity  in  the 
CW  or  quasi-CW  case  and  the  condition  that  the  cross-flow  not  even  halve 
the  gain  coefficient  [see  Equation  (7.2)]  can  be  combined  to  give  the 
condition  that 


gP  <  15  W/cm  , 


(7.9) 


or,  with  the  absorption  coefficient  a  about  3  10'6  cnr*,  that  the  input 
power  satisfy  the  condition 


P  <  5  MW  . 


(7.10) 


The  input  power  for  the  COS  experiments  is  not  to  exceed  5  MW.  The 
range  of  f/number  is  to  be  15  to  30,  which  is  consistent  with  the  condition 
stated  above  for  the  input  power. 

The  CW  or  quasi-CW  case  for  self-focusing  is  not  greatly  affected  by 
subsonic  cross-flow.  Self-focusing  in  xenon  at  HF  wavelengths  is  likely  to 
be  mostly  el ectrostricti ve.  The  gas  pressure  increases  to  balance  the 
el ectrostricti ve  pressure.  The  resulting  increase  in  gas  density  and  in 
refractive  index  is  the  same  as  without  the  cross-flow  to  terms  of  order 
squared  in  the  Mach  number. 

The  pulsed  cases  require  a  different  condition  on  the  input  power. 
Namely,  a  condition  that  limits  the  change  in  the  intensity  distribution 
during  the  SBS  acoustic  wave  rise  time.  Otherwise  the  SBS  phase  conjugation 
would  not  be  able  to  keep  up  with  the  beam  changes.  A  simple  condition 
taken  here  is  that  the  change  in  the  distortion  parameter  be  small  in  a 
phonon  rise  time. 

This  condition  is  written  as 


r  dDP  <  ! 

1  Ht  \  J  < 


(7.11) 


where  r  is  the  phonon  rise  time  and  d(DP)/dt  is  the  time-derivative  of  the 
distortion  parameter  DP.  Upon  evaluating  this  expression  for  a  Gaussian 
beam,  the  condition  becomes 


dA  Tt_  J3_  , 

dt  ,  2  '  1 


(7.12) 


fcV  S'*  (***  *  *  M  *  »  *  >  *  ■ 


or 


60  W/cm 
(t/ts)2 


(7.13) 


where 

r  =  phonon  rise  time  =  0.26  p s  for  40-atm  Xe 
F /#  =  focal  length/beam  diameter 
C$  =  sound  speed 

X  =  wavelength 


(7.14) 


P  <  20  MW/ (t/t 


r,  when  a  =  3  x  10  cm"1 


(7.15) 


The  square  of  the  ratio  of  pulse  time  to  sonic  transit  time  appears, 
as  discussed  earlier.  This  ratio  is  supposed  to  be  less  than  unity.  A 
rectangular  pulse  in  time  is  assumed  for  the  sake  of  simplicity. 

Similar  analysis  is  done  with  pulsed  self-focusing.  The  product, 
phonon  rise  time  times  the  time  derivative  of  the  distortion  parameter 
(actually,  its  negative),  is  to  be  less  than  unity.  The  condition  on  input 
power  is  as  follows: 


f  1  , 

c  V 


Pc  =  ncX  / (128  n2)  ~  1  MW  , 


(7.16) 


where  the  nonlinear  refractive  index  coefficient  n^  is  taken  as 


'  P  ~ 1 

U  =  - d-l  —  :  0.15  x  10  ^  esu,  for  40-atm  Xe, 

16*  p  n 


(7.17) 


e  -  dielectric  constant 
p  -  density 


The  condition  on  input  power  for  the  pulsed  self-focusing  case  is  also 
consistent  with  the  upper  limit  of  5  MW  and  f/number  range  of  15  to  30. 

Optical  breakdown  in  the  SBS  cell  might  be  most  influential.  Experi¬ 
ments  at  1.35  pm  wavelength  (Ref.  7.3)  indicate  a  breakdown  threshold 
intensity  of  about  20  GW'cm2  at  40  atm  xenon.  Scaling  to  HF  wavelength 
(2.8  /im)  by  an  inverse  or  an  inverse-square  wavelength  dependence  yields  an 
estimate  of  5  to  10  GW/cm2  for  the  breakdown  threshold  intensity.  Then, 
the  condition  that  the  pump  intensity  that  is  transmitted  through  the 
forward-posi tioned  SBS  region  to  the  focal  region  be  less  than  this 
threshold  is 

P  <  4  (F/#)2/(i_r) ,  in  kW,  (7.19) 

where  R  is  the  SBS  reflection  coefficient.  The  full  state  of  affairs  for 
the  SBS/breakdown  interaction  is  not  yet  known. 

The  opposing  effects  of  thermal  blooming  and  self-focusing  can  be  made 
to  cancel  with  the  appropriate  absorption  coefficient  and  corresponding 
exponentially  rising  time  history  for  the  pulsed  case  or  exponentially  grow 
in  spatial  profile  in  the  flow  direction  in  the  region  just  forward  of 
the  Rayleigh  range.  The  1/e-folding  time  is  given  by 

P  ^)/[a(7-Dc]  (7.20) 

c  =  1 ight-speed 
7  =  ratio  of  specific  heats 


7.3  M. A.  Greimer-Mothes  and  K.J.  Witte,  preprint,  1986. 
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and  the  1/e-folding  distance  is  the  cross-flow  speed  times  this  time- 
constant  (Ref.  7.4). 

For  xenon  at  40  atm  with  less  than  0.1  ppmv  impurity  level,  the  time- 
constant  is  a  convenient  1  ps.  With  20  m/s  flow  speed,  the  distance  is 
20  /*m,  which  is  to  be  compared  to  the  focal  spot  diameter  of  6  F/*  pm  at  HF 
wavelength. 

7.3  SUBSCALE  EXPERIMENTS 

The  overall  objective  of  the  subscale  experiments  task  has  been  to 
provide  preliminary  experimental  data  to  support  the  PALS  design  activities 
on  phase  conjugation  related  issues.  Specifically,  these  experiments  were 
designed  to  study  the  SBS  threshold  and  reflectivity,  and  conjugation 
fidelity  as  a  function  of  the  following  experimental  parameters:  (1)  f/ 
number  of  the  SBS  telescope,  and  (2)  bandwidth. 

The  f/number  of  the  SBS  telescope  has  a  significant  impact  on  the 
APACHE  system  design.  Based  on  modeling  and  numerical  analysis  it  is 
believed  that  the  threshold  for  the  SBS  process  should  be  independent  of 
the  f/number  of  the  focusing  telescope,  but  prior  to  these  experiments  no 
documented  experimental  results  were  available.  Furthermore,  the  effect  of 
f/number  on  conjugation  fidelity  was  an  issue  which  had  received  relatively 
little  attention  but  significantly  impacts  system  design. 

The  bandwidth  of  the  pump  laser  also  impacts  the  conjugation  process. 
The  SBS  process  requires  a  finite  length  over  which  the  Stokes  beam 
interacts  with  the  pump  to  produce  the  phase  conjugate  return.  If  the 
coherence  length  of  the  laser  is  less  than  this  SBS  interaction  length  then 
the  reflectivity  and  fidelity  of  the  conjugate  return  will  be  degraded.  In 
the  subscale  experiments,  the  plan  was  to  run  the  laser  both  narrow  band 
(with  the  SAM  single  axial  mode  etalon  in  place)  and  broadband  (with  the 
SAM  removed)  with  variations  in  the  SBS  interaction  length  being  controlled 
by  the  f/number  of  the  focusing  lens. 


7.4  P.R.  Longaker  and  M.M.  Litvak,  Journ.  Appl .  Phys.,  40,  4033-4041 
(1969),  see  eq.  (A6) .  Note  on  Eq.  (A3):  Divide  r.h.s.  by  (7-I). 
Note  on  eq.  (A5) :  c  is  light  speed  here  only. 


A  key  output  of  this  experimental  program  is  a  data  base  for  vali¬ 
dating  the  upgraded  capabilities  of  the  BOUNCE  code,  which  was  modified 
during  the  PALS  program  to  include  the  effects  of  near-field  aberrations 
and  absorption.  The  data  also  were  to  supply  information  for  setting 
medium  purity  requirements  on  the  conjugating  material  as  well  as  flow 
rates  in  the  flowing  SBS  cell,  and  to  assist  in  setting  system  level 
requirements  (notably  the  f /number  of  the  SBS  telescope). 

Figure  7-2  shows  a  schematic  layout  of  the  test  setup  used  for  these 
experiments.  A  photograph  of  the  facility  is  presented  in  Figure  7-3.  The 
laser  source  is  a  Molectron  Nd:YAG  laser  operating  at  1.06  with  a  pulse 
length  of  approximately  15  ns  (FWHM)  and  a  maximum  pulse  energy  of  250  mJ 
(in  practice,  to  avoid  damage  to  the  spatial  filter  the  pulse  energy  was 
kept  below  about  100  mJ).  The  beam  hits  a  turning  flat  and  is  directed  to 
a  multielement  telescope  with  a  3  m  effective  focal  length  with  a  spatial 
filter  positioned  at  focus.  This  focal  length  was  chosen  to  keep  the 
fluence  levels  at  focus  below  damage  levels  of  the  spatial  filter  and  air 
breakdown  levels.  Beam  splitters  are  located  on  both  sides  of  the  spatial 
filter  for  joule-meter  readings.  Beyond  the  telescope,  the  beam  is 
recollimated  and  directed  to  the  SBS  focusing  lens  which  is  varied  over  the 
f/number  range  of  10  to  50.  Two  SBS  cells  were  used  in  the  study  using 
CC 1 4  and  xenon.  All  of  the  successful  experiments  performed  to  date 
utilized  the  liquid  CC 1 4  cell.  Further  work  using  the  high  pressure  xenon 
cell  is  to  be  deferred  to  the  APACHE  program. 

Measurements  of  the  characteristics  of  the  laser  are  presented  in 
Figures  7-4  and  7-5.  The  far-field  intensity  distribution  of  the  laser 
measured  using  a  CID  camera  and  the  beam  quality  are  shown  in  Figure  7-4. 
Using  the  power-through-the-pinhole  technique,  the  beam  quality  of  the  beam 
incident  on  the  SBS  cell  is  measured  to  be  BQ= 1.10.  Figure  7-5  shows  the 
bandwidth  characteristics  of  the  laser  measured  using  a  Fabry-Perot  inter¬ 
ferometer.  With  the  SAM  (single  axial  mode)  etalon  in  place,  the  bandwidth 
of  the  laser  at  1.06  pm  is  measured  to  be  approximately  200  MHz,  yielding  a 
coherence  length  of  1.5  m.  The  P(t)  trace  obtained  simultaneously  with  the 
bandwidth  measurements  shows  the  characteristic  shape  of  a  Q-switched  laser 
with  a  sharp  rise  followed  by  a  somewhat  longer  decay.  With  the  SAM 
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CIO  Isometric 


Contour  Plot 


BQ  into  SBS  cell  =  1.10 

•  Power  through  pinhole  method 

Figure  7-4.  Measurements  of  the  Nd:YAG  laser  far-field  intensity 
distribution  and  beam  quality. 

removed,  the  character  of  the  P(t)  trace  is  significantly  altered,  with 
several  peaks  and  dips  noted  in  the  trace.  The  Fabry-Perot  inteferogram 
for  this  case  shows  the  bandwidth  to  be  about  13  GHz,  corresponding  to  a 
coherence  length  of  slightly  over  2  cm. 

Several  examples  of  SBS  conjugation  using  P(t)  traces  from  the  CCI4 
liquid  cell  are  presented  in  Figure  7-6.  In  the  oscilliscope  traces,  both 
the  incident  pump  beam  and  the  Stokes  return  (time  delayed  by  approximately 
50  ns  on  the  trace)  are  shown  on  the  same  sweep.  The  first  series  of  tests 
shows  the  effect  of  amplitude  of  the  pump  beam  on  the  conjugate  return. 

For  an  incident  pump  beam  energy  of  1.9  mJ,  no  Stokes  return  is  observed. 
This  result  is  consistent  with  the  calculated  energy  threshold  for  SBS  in 
CCI4  of  4.4  mJ  (based  on  a  conversion  from  power  to  energy  using  FWHM) . 

When  the  pump  energy  is  increased  to  2.9  mJ,  a  Stokes  return  is  observed 
with  a  peak  reflectance  of  approximately  50  percent,  while  reflectivity 
measured  in  terms  of  energy  is  14  percent.  This  disparity  between  power 
and  energy  reflectivity  results  from  the  fact  that  the  phonon  rise  time  in 
CCI4  at  1.06  /*m  (9  ns)  is  comparable  to  the  duration  of  the  pump  beam 
pulse.  Note  that  the  observance  of  an  SBS  return  at  an  energy  level  below 
the  calculated  value  is  an  indication  that  the  computation  of  threshold 


energy  using  the  pulse  FWHM  to  convert  from  power  to  energy  is  only  an 
approximation .  A  more  definitive  comparison  would  be  to  compare  the 
instantaneous  power  form  the  P(t)  trace  at  the  point  in  time  when  the  SBS 
process  begins,  and  such  a  determination  is  presently  being  made.  In  the 
last  trace  of  this  series,  the  pump  beam  energy  is  increased  to  17.1  mJ, 
and  the  reflectivity  of  the  Stokes  return  is  about  50  percent  in  terms  of 
energy  and  over  80  percent  in  terms  of  peak  power. 

Also  shown  in  Figure  7-6  is  an  example  of  a  trace  obtained  when  the 
laser  was  operating  multimode,  such  that  evidence  of  mode  beating  is 
observed.  The  interesting  feature  from  these  results  is  that  the  high 
frequency  structure  of  the  mode  beating  which  is  superimposed  on  the  pulse 
is  reproduced  by  the  conjugate  return,  even  though  the  rise  time  for  this 
structure  (approximated  to  be  less  than  1  ns)  is  significantly  faster  than 
the  phonon  rise  time.  This  type  of  response  is  predicted  for  broadband 
input  as  long  as  the  coherence  length  is  longer  than  the  interaction 
length.  Some  attenuation  in  this  high  frequency  structure  is  observed  in 
the  SBS  return,  however,  indicating  that  the  phonon  rise  time  has  acted  as 
a  1 ow  pass  f i 1  ter. 

The  remaining  experimental  data  to  be  presented  deal  with  the  effects 
of  f/number  on  conjugation  threshold  and  reflectivity.  Figure  7-7  shows 
the  results  of  SBS  energy  reflectivity  experiments  performed  as  a  function 
of  f/number  of  the  SBS  focusing  lens  over  the  range  of  f/10  to  f/50.  The 
data  are  plotted  in  the  form  of  reflected  energy  as  a  function  of  incident 
beam  energy.  From  these  results  it  is  observed  that  both  the  threshold  and 
reflectivity  of  the  SBS  conjugate  return  (having  a  slope  efficiency  of  90 
percent)  are  independent  of  f/number  over  the  range  of  f/numbers 
investigated. 

An  interpretation  of  these  results  may  be  made  using  the  expression 
for  the  small  signal  gain.  From  analytical  modeling  it  is  known  that  the 
threshold  for  SBS  phase  conjugation  occurs  when  the  SBS  small  signal  gain, 
equal  to  g,  the  gain  coefficient  of  the  medium,  times  I,  the  intensity, 
times  1,  the  pump'Stokes  interaction  length,  is  greater  than  approximately 
25.  As  the  f/number  of  the  focusing  lens  changes,  both  the  intensity  of 
the  beam  near  focus  and  the  depth  of  focus  (OOP)  change  such  that  I (DOF)  is 


Figure  7-6.  Examples  of  SBS  reflectivity  using  CCI4. 
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SBS  threshold  and  reflectivity  independent  of  f/number 

•  Slope  efficiency  =  0.9 

•  Threshold:  3.5  mJ  (measured) 

4.4  mJ  (calculated  based  on  FWHM) 

Figure  7-7.  Effect  of  f/number  on  reflectivity 

a  constant.  The  data  shown  in  Figure  7-7  suggest  that  as  the  intensity 
varies,  the  length  of  the  region  over  which  significant  interaction  between 
the  pump  and  the  Stokes  correspondingly  changes  such  that  the  gIL  product 
exceeds  that  for  conjugation  threshold.  It  is  believed  that  this  process 
will  continue  until  the  required  interaction  length  for  SBS  exceeds  the 
coherence  length  of  the  laser  after  which  both  the  reflectivity  and  the 
fidelity  of  the  conjugate  beam  will  be  degraded. 

The  data  presented  in  Figures  7-7  summarize  the  present  status  of  the 
experiments  performed  on  the  PALS  program.  As  part  of  the  follow-on  work 
to  be  performed  on  the  APACHE  program,  several  additional  experiments  are 
planned  to  address  all  of  the  issues  previously  identified  in  the  experi¬ 
mental  objectives.  This  work  will  include  examining  the  effects  of 
f-number  on  conjugation  fidelity,  aberrations  in  the  near  field,  absorption 
in  the  conjugating  medium,  and  effects  of  finite  phonon  rise  time  using 
xenon  as  the  conjugation  medium. 
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One  of  the  issues  which  still  remained  at  the  conclusion  of  the  PALS 
program  in  the  area  of  subscale  experiments  was  conjugation  using  xenon  gas. 
Preliminary  attempts  to  use  a  xenon  cell  for  phase  conjugation  were  unsuc¬ 
cessful,  in  that  poor  reflectivity  was  obtained,  and  gas  breakdown  may  also 
have  been  occurring.  Following  the  CDR  presentation  at  NRL,  additional  ef¬ 
fort  was  undertaken  to  assess  the  cause  of  these  difficulties  as  part  of  the 
final  report  preparation.  This  work  revealed  that  a  potential  cause  for 
these  problems  lies  in  contamination  of  the  gas  caused  by  impurities  within 
the  cell.  Because  of  the  specific  nature  of  the  design  of  the  cell,  attempts 
to  remove  these  impurities  proved  to  be  unsuccessful.  Since  these  impurities 
appeared  to  be  the  cause  of  the  poor  conjugation  performance,  a  new  cell  was 
designed  with  medium  purity  as  an  important  requirement.  The  cell  has  been 
fabricated,  and  preliminary  inspection  suggests  that  the  medium  purity  dif¬ 
ficulties  experienced  earlier  will  not  be  encountered  in  the  new  design. 


